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Abstract 

Stochastic calculus with respect to fractional Brownian motion (fBm) has attracted a lot of interest 
in recent years, motivated in particular by applications in finance and Internet traffic modeling. 
Multifractional Brownian motion (mBm) is a Gaussian extension of fBm that allows to control the 
pointwise regularity of the paths of the process and to decouple it from its long range dependence 
properties. This generalization is obtained by replacing the constant Hurst parameter H of fBm by 
a function h{t). Multifractional Brownian motion has proved useful in many applications, including 
the ones just mentioned. 

In this work we extend to mBm the construction of a stochastic integral with respect to fBm. 
This stochastic integral is based on white noise theory, as originally proposed in [15], [6], [4] and 
in [5]. In that view, a multifractional white noise is defined, which allows to integrate with respect 
to mBm a large class of stochastic processes using Wick products. Ito formulas (both for tempered 
distributions and for functions with sub-exponential growth) are given, along with a Tanaka Formula. 
The cases of two specific functions h which give notable Ito formulas are presented. 
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1 Background and Motivations 

Fractional Brownian motion (fBm) [3T] is a centered Gaussian process with features that make it 
a useful model in various applications such as financial and teletraffic modeling, image analysis and 
synthesis, geophysics and more. These features include self-similarity, long range dependence and the 
ability to match any prescribed constant local regularity. Fractional Brownian motion depends on a 
parameter, usually denoted H and called the Hurst exponent, that belongs to (0,1). Its covariance 
function Rh reads: 

i?H(t,s):=^(|f^ + kr-|t-sr), (1.1) 
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where is a positive constant. A normahzed fBm is one for which = 1. Obviously, when H = ^, 
fBm reduces to standard Brownian motion. 

A useful representation of fBm B^^^ of exponent H is the so-called liarmonizable one: 

Cff Jr \u\ ' 

2cob(7ra:)r(2-2x) j _ / ^ — ^ j for X in (0, 1) and W denotes the complex- valued 

x{\ — 2x) j \ i (2a:-t-l) sin(7ra;) / \ ^ I Jr- 

Gaussian measure which can be associated in a unique way to VF, an independently scattered standard 
Gaussian measure on M (see (SS] p. 203-204 and [35] p. 325-326 for more information on the meaning of 
J^f{u) W{du) for a complex-valued function /). From and Gaussianity, it is not hard to prove 

that fBm is i7-self-similar. 

The fact that most of the properties of fBm are governed by the single number H restricts its application 
in some situations. Let us give two examples. The long term correlations of the increments of fBm decay 
as k'-^^~'^\ where k is the lag, resulting in long range dependence when H > 1/2 and anti-persistent 
behavior when H < 1/2. Also, almost surely, for each t, its pointwise Holder exponent is equal to 
H. Since H rules both ends of the Fourier spectrum, i.e. the high frequencies related to the Holder 
regularity and the low frequencies related to the long term dependence structure, it is not possible to have 
at the same time e.g. a very irregular local behavior (implying H close to 0) and long range dependence 
(implying H > 1/2). As a consequence, fBm is not adapted to model phenomena which display both 
these features, such as Internet traffic or certain highly textured images with strong global organization. 
Another example is in the field of image synthesis: fBm has frequently been used for generating artificial 
mountains. Such a modeling assumes that the regularity of the mountain is everywhere the same. This 
is not realistic, since it does not take into account erosion or other meteorological phenomena which 
smooth some parts of mountains more than others. 

Multifractional Brownian motion (mBm) [33l|3] was introduced to overcome these limitations. The basic 
idea is to replace in (|1.2p the real by a function h{t). More precisely, we will use the following definition 
of mBm: 

Definition 1.1 (Multifractional Brownian motion). Let h : R ^ (0,1) be a continuous function and 
a : (0, 1) — > M 6e a function. A multifractional Brownian motion with functional parameters h and 
a is defined as: 

B^^^^'\t) = a{h{t)) [ Widu). (1.3) 

Jr \u\ ^ ' ' 



1 



2 



s 



(1.4) 



Its covariance function reads [5]: 

R{h,a){t,s) ^ a{h{t)) a{h{s)) 

where ht^s '■= and Cx has been defined in (|1.2p . 

It is easy to check that mBm is a zero mean Gaussian process, the increments of which are in general 
neither independent nor stationary. 

For T in R^j,, we will again call (/i, Q;)-multifractional Brownian motion on [0,r] the centered Gaussian 
process whose covariance function is equal to Rf^h.a) on [0,T] x [0,T]. 
When a ~ ac : X 1-^ — , wc get that: 



C{h(t))C{h(s)) 



(1.5) 



As a consequence, if h is the constant function equal to the real H, then is a normalized fBm. For 

this reason, we will call B^^'""^ a normalized mBm. Since in the sequel we will consider only normalized 
mBm, we simplify the notation and write from now on B^^'^ for B^'^'""'^ and Rh for R(h,aa)- 
One can show |18[[T^ that the pointwise Holder exponent at any point t oi B^^^ is almost surely equal to 
h{t) A Ph{t), where I3h{t) is the pointwise Holder exponent of h at t. In addition, the increments of mBm 
display long range dependence for all non-constant h{t) (the notion of long range dependence must be 
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re-defined carefully for non-stationary increments, see Finally, at least when h is C^, for all u G R, 
mBm locally "looks like" fBm with exponent h{u) in the neighbourhood of u in the following sense |33j : 

where the convergence holds in law. These properties show that mBm is a more versatile model that 
fBm: in particular, it is able to mimic in a more faithful way local properties of financial records, Internet 
traffic and natural landscapes [3 1301 114j by matching their local regularity. This is important e.g. for 
purposes of detection or real-time control. The price to pay is of course that one has to deal with the 
added complexity brought by having a functional parameter instead of a single number. 
Because of applications, in particular in finance and telecommunications, it has been an important 
objective in recent years to define a stochastic calculus with respect to fBm. This was not a trivial 
matter, as fBm is not a semi- martingale for H ^ ^. Several approaches have been proposed, based 
mainly on Malliavin calculus [TH[TJ[T7], pathwise approaches and rough paths ([SnUHllIS] and references 
therein), and white noise theory |15[ [51[5]. Since mBm seems to be a more flexible, albeit more complex, 
model than fBm, it seems desirable to extend the stochastic calculus defined for fBm to it. This is the 
aim of the current work. In that view, we will use a white noise approach, as it offers several advantages 
in our frame. The main task is to define a multifractional white noise as a Hida stochastic distribution, 
which generalizes the fractional white noise of, e.g., [151 16]. For that purpose, we use the properties 
of the Gaussian field (B^^^t)) //)gex(o i)- particular, it is a crucial fact for us that the function 
H is almost surely C°°. This entails that multifractional white noise behaves essentially as 

fractional white noise, plus a smooth term. We obtain an Ito formula that reads: 

/(T,B"'>(T)) = /(0,O)+ / ?l{t.B>'->[i))dt+ f '^^(t,B<">{mdB"'>(t) 
where the meaning of the different terms will be explained below. 

The remaining of this paper is organized as follows. In section [51 we recall basic facts about white noise 
theory. We study a family of operators, noted {Mh) fj£[Q i)i which are instrumental for constructing the 
stochastic integral with respect to mBm in section [3] Section |4] defines the Wiener integral with respect 
to mBm. We build up a stochastic integral with respect to mBm in section O Various instances of Ito 
formula are proved in section [HI Finally, section [7] provides a Tanaka formula, along with the study of 
two particular h functions that give notable results. Readers familiar with white noise theory may skip 
the next section. 

2 White noise theory 

We recall in this section the standard set-up for classical white-noise theory. We refer e.g. to [371 [H] for 
more details. 

2.1 White noise measure 

Let := {/ e C°°(M) : V(p, (?) G N^, lim \xP f^'^'>{x)\ = 0} be the Schwartz space. A family of 

|a;|— f-t-oo 

functions {fn)nen (^(K))^ is said to converge to as n tends to +cx) if for all (p, q) in W' we have 
lim sup \x'P fn\x)\ = 0. The topology hence given on y{R) is called the usual topology. Let ^'(M) 

n->-|-oo 

denote the space of tempered distributions, which is the dual space of S^{U). The Fourier transform of 
a function / which belongs to ^-'^(R) U L^(M) will be denoted / or T{f): 

Him m / e-^-^f{x)dx, ^ G M. (2.1) 
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Define the probability space as ft :^ (M) and let := B{S^ (M)) be the tr-algebra of Borel sets. The 
Bochner-Minlos theorem ensures that there exists a unique probability measure on O, denoted ^, such 
that: 

/ e'<"'^>/i(dw) =e"*"^"''(«), V/eJ5^(M), (2.2) 

where < w,/ > is by definition i.e the action of the distribution a; on the function /. For / in 

■y{M.) the map. noted < .,/ >, from fi to M defined by < .,/ > (i^) =< uj,f > is thus a centered 
Gaussian random variable with variance equal to ||/||^2(jj-| under the probability measure fj,, which is 
called the white-noise probability measure. In other words, 



IE[<.,/>]=0, E[< .,/>2] = 11/ 



for all / in S^{M.). Besides, for a measurable function F, from ^ (R) to R, the expectation of F with 
respect to fi is defined, when it exists, by E[F] := Ef^[F] := J^^ F {w) ^(dui) . Equality p.2p entails that 
the map C defined on ,_5^(R) by 

C: (^(R),<,>i2(R)) ^ {LH^,^,t^),<,>{mn,^,^.)) on 
/ ^ C(/) :-<.,/> 

is an isometry. Thus, it extends to L^(M) and we still note C this extension. For an arbitrary / in L^(R), 
we then have < .,/ >:= lim <.,/«> where the convergence takes place in L'^{i},J^, fj.) and where 

n— ^+00 

(/n)„gN ^ sequence of functions which belongs to o5^(M) and converges to / in i^(R). In particular, 
define for all t in R, the indicator function il[o,t] by 

1 i/ 0^ss:t, 
il[o.t](s) = { ^ if t ^ s ^ except if < = s = 

otherwise, 

Then the process (-Bt)^^^, defined for t e R, on f2 by: 

Bt{ij) ■.^B{t,uj) := <Lj,il[o^t]> (2.4) 

is a standard Brownian motion with respect to fi. It then admits a continuous version which will be 
denoted B. Thanks to (|2.4p we see that, for all functions / in L^(R), 

Iiif){Lu) = <LuJ>= f f{s)dBs{uj) 11 - a.s. (2.5) 

2.2 Properties of Hermite functions and space =5^'(M) 
For every n in N, define 

2 ' 

/i„(x) := (-l)"e'^ — (e"^') the n*'* Hermite polynomial (2.6) 

and 

en{x) ■.= TT-^''^{2''n\)~^''^e-''^''^K{x) the n*'' Hermite function (2.7) 
We will need the following properties of the Hermite functions: 

Theorem 2.1. 1. The family (efc)j,gpj belongs to and forms an orthonormal basis of L'^{R) 

endowed with its usual inner product. 

2. There exists a real constant C such that, for every k in N, maxje^ (x)| C (fc + l)"^/^^ More 
precisely, there exist positive constants C and 7 such that, for every k in N, 

' ^ ^' \Ce'-"= if \x\ > 2VkTT. 
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See [37| for proofs. 

In order to study precisely J^(E) and its dual ,5^'{K) it is desirable to have a family of norms on the 
space o5^(R) which gives us the usual topology. 

Definition 2.1. Let (| \^) ^ he the family oj norm defined by 



+00 

|/|^:=5](2fc + 2f^</,efe>i.(K), V(p, /) G Z x l2(M). (2.9) 

fc=0 

For p in N, define the spaces ^p(R) :~ {/ G i^(M), < +00} anrf J?^_p(]R) as 6em^ t/ie completion 
of the space (M) with respect to the norm \ | _^ . 

Remark 2.2. For a function f which is not in L^(M), we may still define \f\^ by allowing \ f\^ to he 
infinite. 

It is well known (see [22]) that the Schwartz space ,5^{R) is the projective limit of the sequence {S^p{M))^^^ 
and that the space ,5^' {M.) of tempered distributions is the inductive limit of the sequence {5f.p{E^) p^-^. 
This means firstly that both equalities = n ^p(M) and ^'(R) = U ^p(R) hold. Secondly that 

convergence in J5^(M) is nothing but convergence in J?'p(R) for every p in N and that convergence in 
is convergence in ^p(R) for some p in N. Moreover one can show that, for any p in N, the dual 
space oS^p(R) of o5^p(R) is ^p(R). This is the reason why we will write ^p(R) in the sequel to denote the 
space ^p(R). Finally one can show that the usual topology of the space =5^(R) and the topology given by 
the family of norms (| j^) arc the same (see [20j appendix A. 3 for example). Moreover, convergence 

in the inductive limit topology coincides with both convergence in the strong and the weak* topologies 
of J5^'(R). 

In view of definition 12. 11 it is convenient to have an operator defined on =5^(R) whose eigenfunctions are 
the sequence (e„)^gpj and eigenvalues are the sequence (2n + 2)„gN. It is easy to check that the operator 
A, defined on J5^(R), by 

A:=-^+x2 + l (2.10) 

verifies these conditions. Moreover A is an isometry from (.^^(R), | |-^) to (L^(R), | |g). 

It may thus be extended to an operator, still denoted A, from o$^(R)' = =5^i(R) to L^(R). It is then easy 
to show (see [57] p. 17-18 for instance) that A is invertible and that its inverse A~^ is a bounded operator 
on (R) . Let us note | g | q : = 1 1 <? 1 1 ^2 (rj for any g in (R) . For p in Z, let A^ denote the p*^ iteration of the 
operator A, if p belongs to N, and of A~^ otherwise. Then Bom{AP) ,yp{R) and Bom{A-P) L'^{R) 
where ]D)om(t/) denotes the domain of the operator U and p belongs to N. Moreover, for every g in Z and 
every / := J2t=o < f^(^k >mR)ek in Dom(A9), the equality A^f = 2+^^ (2fc + 2)' < /, >mR)ek 
holds. Hence, 

+00 

l/l' = I^VIo = E(2^ + 2)2'</'^^->i=(R)' ^9^^- (2.11) 

k=0 



2.3 Space of Hida distributions 

From now on we will denote as is customary (L^) the space L'^{U,0, p) where G is the cr-field generated 
by (< -J / >)/gl2(r)- Neither Brownian motion nor fractional Brownian motion, whatever the value of H , 
are diff'erentiable (see [31] for a proof). However, it occurs that the mapping t H> B'^^\t) is differentiable 
from R into a space, noted (5)*, called the space of Hida distributions, which contains [L"^)- In this 
section we recall the construction of {S)* . 

For every random variable <& of (L^) there exists, according to the Wiener-Ito theorem, a unique sequence 
(/n)„gN of functions /„ in _L^(R") such that $ can be decomposed as 

+00 

* = E^"(/») (2.12) 

n=0 
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where L^(]R") denotes the set of all symmetric functions / in L^(R") and denotes the n^^ multiple 

Wiener-Ito integral of / defined by 

/„(/):= / f{t)dB'\t)^n\f{f"---{f'f{h,---,t^)dB{h))dB{t2)---dB{tn)), (2.13) 

with the convention that /o(/o) = /o for constants /q. Moreover we have the isomctry = 
+00 . 

2 « 2 

^^n!||/„|| 2,2(R>.)- For any $ := ^ In{.fn) satisfying the condition ^ n\ |A'^"/„|q < +00, define 

r. n— n— 

n— 

the element r(A)($) of (L^) by 

r(A)($):=^/„(^«"/„), (2.14) 

n=0 

where A*^" denotes the n*'' tensor power of the operator A defined in (|2.10p (see [IS] appendix E for 
more details about tensor products of operators). The operator V{A) is densely defined on (L^) and 
is called the second quantization operator of A. It shares a lot of properties with the operator A. In 
particular it is invertible and its inverse T{A)^^ is bounded (see [23). Let us denote \\^p\^ := || 95 11(^2) 
for any random variable <p in (L^) and. for n in N, let Dom(r(A)") be the domain of the n"^ iteration 
of T{A). The space of Hida distributions is defined in a way analogous to the one that allowed to define 
the space ^'(R): 

Definition 2.2. Define the family of norms (|| || ) by: 

:= l|r(A)f$||o = ||r(A)p$||(^.), Vpez, v$ e (l^) n©om(r(A)P). (2.15) 

For any p in N, let (Sp) := {$ € (i^) : T{A)p^ exists and belongs to (i^)} and define (5_p) as being 
the completion of the space (i^) with respect to the norm \ \ \\_p. 

As in [27], we let (S) denote the projective limit of the sequence {{Sp))^^^^ and (S)* the inductive 
limit of the sequence {{S-p))^^^^. Again this means that we have the equalities (S) ~ n^{Sp) (resp. 

(S)* = U (5_p)) and that convergence in (S) (resp. in (S)*) means convergence in (Sp)* for every p 

in N (resp. convergence in (S-p) for some p in N ). The space (S) is called the space of stochastic test 
functions and (S)* the space of Hida distributions. As previously one can show that, for any p in N, 
the dual space (Sp)* of Sp is (S-p). Thus we will write (S-p), in the sequel, to denote the space (Sp)* . 
Note also that (S)* is the dual space of (S). Wc will note < , >■ the duality bracket between (S)* and 
(5). If $ belongs to (L^) then we have the equality <:^,ip:> = < ^,(p >(l^) = Furthermore, 
as one can check, the family (|/|p) is an increasing sequence for every / in J5^(M). Thus the family 

(II < ., / > II ) is an increasing sequence for every / in y{R). 

Remark 2.3. (i) A consequence of the previous subsection is that for every element f :— anCn 
in oS^'(K) where (a„)^gf^ belongs to R^, there exists po in N such that f belongs to c5fp„(M). Moreover 
if we define <f> := '^n=a < •,e„ >, then $ belongs to (S-pg) C (S)* and we have W^lfLpg ~ l/l-po ~ 
J2n^o (27t+2)^p" +0O- Conversely, every element written as $ X]n=o < ■,e„ > where {bn)„^ff 
belongs to , belongs to (S)* if and only if there exists an integer po such that "^^S) (2rt+2)^''" ^ +00. 
In this case the element f := "^^=0 bn^-n belongs to ^po(R) and hence is a tempered distribution which 
verifies \ftp„ = mtp„ = En=o (2„+"2)^''° ■ 

(m) Let ^(R") be the space of symmetric Schwartz functions defined on R". Let p in N* and $ be 
an element of (S-p). Then 4> can be written $ := J^nS) In{Fn)t where {Fn)^^^ is a sequence in 

(^'(R"))'^, where i^'(R") is the dual o/i^R"). 

Since we have defined a topology given by a family of norms on the space [S)* it is possible to define a 
derivative and an integral in iS)* (see [H] chapter 3 for more details about these notions). Let / be an 
interval of R (which may be equal to R). 
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Definition 2.3 (stochastic distribution process). A function $ : / — > (S)* is called a stochastic 
distribution process, or an [S]*^ process, or a Hida process. 

Definition 2.4 (derivative in {S)*). Let to E I. A stochastic distribution process $ : / — > (S)* is said 
to be differentiate at <o */ the quantity lini r^^ ($(^0 + r) — *i'(^o)) exists in {S)* . We note ^(to) 

the {S) -derivative at to of the stochastic distribution process $. <I> is said to be differentiable over I if it 
is differentiable at to for every to in I . 

The process $ is said to be continuous, C^, ■ ■ ■ ,C^, ■■ ■ in [S)* if the (iS)*-valued function $ is, con- 
tinuous, C^, • • • , C'^, • • • . We also say that the stochastic distribution process $ is (5)* -continuous and 
so on. It is also possible to define an (5)* -valued integral in the following way ([27l [22]). We first 
recall that (R, dt) denotes the set of measurable complex- valued functions defined on K. such that 



Theorem-Definition 2.1 (integral in (S)*). Assume that $ : R — > (S)* is weakly in L^{M.,dt), i.e 
assume that for all (p in (S), the mapping u i— ><C ^{u), tf~> from R to R belongs to L^{M.,dt). Then 
there exists an unique element in (S)* , noted J„^(u)du such that 



^{u)du,Lp'> = J <$(u),v3> du for all if in (S) . (2-16) 

We say in this case that $ is (iS)*-integrable on R in the Pettis sense. In the sequel, when we do not 
specify a name for the integral of an (5)*-integrable process $ on R, we always refer to the integral of 
$ in Pettis' sense. See [17] p.245-246 or [H] def. 3.7.1 p.77 for more details. 

Recah from 1^ that /i(/) = /(s) dB{s). For p in N, let I'f^ : (^^(R), | \^^) {{L^),\\ ||_p) be 
the map defined by l[''\f ) ~< ., / >. For / in L^{R) and p in N we have, using [22 p 26, the equality 

wii'Hm-, = mAr^ihifm, = mA-pifm, = \A-pif)\, = i/u^. 

Hence is an isometry and we can extend it to ^_p(R). Since we can do this for every integer p in N, 
we may then give a meaning to Ii{f) ~< ., / > for every tempered distribution / as being the element 
4^^(/) where / belongs to J5^-p(R). 

Remark 2.4. We will also note L f{s) dB{s) the quantity < .,/ > when f belongs to S"{R). Hence 
we give a meaning to the quantity Jj^ f{s) dB{s) for every f in ^'(R). 

2.4 S-transform and Wick product 

For rj in =5^(R), the Wick exponential of < .,rj >, denoted : e^ '''^ :, is defined as the element of [S) 
given by : e<-''> : =^ Yl=n ^^"^ h{r]®^) (equality in (L^)). More generally, for / e L^{B), we define 
■ e< J> : as the (L^) random variable equal to e^ ''^^~il'^lo (see [5S] theorem 3.33). We will sometimes 
note exp* < ., / > instead of : e"^ '^^ :. This random variable belongs to LP{rt, n) for every integer p ^ 1. 
We now recall the definition of the ^-transform of an element $ of (S*), noted S'($) or S'[$]. 5'($) is 
defined as the function from ^(R) to R given by 

V77e^(R), 5'($)(77) ;= <$,: e<-''> :> . (2.17) 

Note that S^{ri) is nothing but E[$ : e<-'''> :] = e^^l^'lo E[$ e< '''>] when $ belongs to (L^). Following 
[5], formula (6) and (7), define for 77 in ^(R) the probability measure Q,, on the space (0,J^) by its 

Radon-Nikodym derivative given by e^ '''^:. The probability measures Q,, and ^ are equivalent. 

Then, by definition, 

\f^e{L'), 5($)(77)=EqJ$]. (2.18) 
Lemma 2.5. 1. Let p be a positive integer and ^ be an element of (S-p). Then 

^ ll^ll-p e*!''!?, for any rj m J^(R). (2.19) 
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2. Let $ := cik < i^k > belong to (5)*. The following equality holds for every i] in ^(M); 

+ CSO 

WW = 5Z«fc<'7,e;. >i2(R). (2.20) 

fc=0 

Proof. Item 1 is proved in [27] p. 79. Item 2 is an easy calculation left to the reader. □ 
Another useful tool in white noise analysis is the Wick product: 

Theorem-Definition 2.2 ([27J p. 92). For every ($, *) e (5)* X {Sy, there exists a unique element of 
(S)* , called the Wick product of $ and 'i' and noted $ o vj/^ such that, for every rj in ^(R), 

S{^o^){^) = S{^){r,)S{mri)- (2.21) 
Lemma 2.6. For any {p,q) € N'^ , X e (5_p) and Y e (S-q), 

\SiXoYm < \\X\\_^ \\Y\\_^ el^l-^'-t-). (2.22) 
Proof. The proof is easy since, for every 77 in ^(M), using lemma [^31 

\S{XoYM\ = \SiX)i^)\ \SiY)i7^)\ ^ \\X\\_^ eil< ||r||_, e*!"!? < \\X\\_^ \\Y\\_^ gl'^lL.,.;,) . □ 

k times 

For any $ in (S)* and k in N* let ^^'^ denote the element $ o • • • o $ of (S)* . We can generalize the 
definition of exp* to the case where <& belongs to (S)*: 

Definition 2.5. For any $ in (S)* such that the sum ^;|- converges in (5)*, define the element 

expO $ of {Sr by cxp- <D ^+^0 ^ . 

For / in L^(R) and $ :=< .,/ >, it is easy to verify that exp* $ given by definition 12.51 exists and 
coincide with : e^ '^^ : defined at the beginning of this section. 

Remark 2.7. // $ is deterministic then, for all in (S)* , <f> o 5* = ^'I'. Moreover, let (^t)tgjj be a 

Gaussian process and let % he the subspace of (L^) defined by TL := vectB.{Xt]t G M} . If X and Y 
are two elements ofH then X oY ^ XY - E[XY]. 

We refer to [23 chapters 3 and 16 for more details about Wick product. The following results on the 
S-transform will be used in the sequel. See [27j p. 39 and [5T] p. 280-281 for proofs. 

Lemma 2.8. The S-transform verifies the following properties: 

(i) The map S : ^ y-^ from [S)* to [S)* , is infective. 

(ii) Let $ : M — 5- {S) be an [S) process. If ^ is (S) -integrable over M then 

S{ I $(w) du){?]) = / S'($(u))(?7) du , for all rj m ^(R). 

(Hi) Let $ : M — > [S) be an (S) -process differentiable at t. Then, for every 77 in S^(R) the map 
u [S^{u)]{ri) is differentiable at t and verifies 5'[^(t)](ry) = ^ [5'[$(t)](r/)] . 

It is useful to have a criterion for integrability in {S)* in term of the 5-transform. This is the topic of 
the next theorem (theorem 13.5 in j27j). 

Theorem 2.9. Let $ : R — > (S)* be a stochastic distribution process satisfying: 
(i) The map t i— > S[^{t)\{vi), from R to R, is measurable for all rj in ,_5^(R). 

(ii) There is a natural integer p, a real a and a function L in L^{R,dt) such that for all 77 in c5^(R), 

\s{m){v)\^m e'^K. 

Then $ is (S)* -integrable over R. 
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Lastly, when the stochastic distribution process is an (L^)-valued process, the fohowing resuh holds (see 

my- 

Theorem 2.10. Let X : M — > (i^) be such that the function t i— >■ S{Xt){r]) is measurable for all 77 in 
J^{M.) and that 1 1— > ||^t||o is in L^{R,dt). Then X is (S)* -integrable over M and 



[ Xtdt < / ll^tllo^^- 

JR "'R 



3 The operators Mh and their derivatives 

3.1 Study of Mj:^ 

Let us fix some notations. We will still note u or the Fourier transform of a tempered distribution 
u and we let Ljg^{M.) denote the set of measurable functions which are locally integrable on M. We also 
identify, here and in the sequel, any function / of with its associated distribution, also noted Tf. 

We will say that a tempered distribution v is of function type if there exists a locally integrable function 
/ such that V = Tf (in particular, < v,<f) >= f{t) (l){t) dt for in ^(R)). 

Let H € (0,1). Following [T3], we want to define an operator, denoted M//, which is specified in the 
Fourier domain by 



M^){y) ^ \y\"^-"u{y). y G K*. (3.1) 

CH 

This operator is well defined on the homogeneous Sobolev space of order 1/2 — H, L|^(R): 

L%{M.) := {u e J^'(M) ; u = Tf; f e L]^M) and < +00}, (3.2) 

where '■— -p- derives from the inner product on L|^(M), defined by: 



<u,v>h:=^ I \i\''^"u{0v{0dt (3.3) 

and ch has been defined right after formula (|1.4p (the normalization constant will be explained in 
remark [?75|) . It is well known - see [5] p. 13 for example - that (L|j(M), <, >^) is a Hilbert space. The 
nature of the spaces (R) when H spans (0, 1) is described in the following lemma, the proof of which 
can be found in j8] pl5, theorem 1.4.1 and corollary 1.4.1. 

Lemma 3.1. If H is in (0,1/2], the space L|^(R) is continuously embedded in L^^^(R). When H is in 
[1/2,1), the space L^/^(M) is continuously embedded in L|j(R). 

Since Mh{u) belongs to L^(R) for every u in L|j(R), Mh is well defined as its inverse Fourier transform, 
i.e.: 

Mh{u){x) — F{Mh{u)){—x), for almost every a; in R. (3-4) 
27r 

The following proposition is obvious in view of the definition of Mh: 

Proposition 3.2. Mh is an isometry from (L|^(R), <, >^) to (L^(R), <, >^2(][{-|). 

Let £{M.) denote the space of simple functions on R, which is the set of all finite linear combinations of 
functions il[a,&](-) with a and b in R. It is easy to check that both ^(R) and £(R) are subsets of L|f(M). 
It will be useful in the sequel to have an explicit expression for Mnif) when / is in ^(R) or in £(]R). 
To compute this value, one may use the formulas for the Fourier transform of the distributions | |", a in 
(— 1, 1), given for instance in [3] (chapter 1, § 3). This yields, for almost every x in R, 



MHill[aM)i^) - 2cHr{H+l/2) cos(f (if~l/2)) 



H_^|3/2-H |„_3,|3/2-H 



(3.5) 



By the same method, for / in ^(R) one gets, for almost every real x: 
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for < i/ < 1/2 MHif){x) = 7H / ^ ^j^-A^dt (3.6) 



f{x-t)-f{x) 



fori/ = 1/2 = /(x) (3.7) 

fA 

\t-x\ 



for 1/2 < if < 1 MH{f){x) = -iH / —^^^dt (3.8) 



/ . \ i 

where 7^ := 2cHr(H-i/2^s(f (g-1/2)) = 2r(Hi'vi cos'(f"g-i/2)) ■ When / belongs to ^(R), equahty 
(A.l) of [TS] yields (up to a constant), for almost every x in R, 

d 



MH{f){x) ^an- 

dx 



{t~x)\t-xf~^'^f{t)dt 



(3.9) 



"^^"^^ = 2c«r(if + l/2)coI(f(i/-l/2))- ('-^"^ 

In order to extend the Wiener integral with respect to fBm to an integral with respect to mBm (in 
section |4?2|) we will need the following equality: 



Proposition 3.3. £{M) ' " = L']j{M). 

This is a straightforward consequence of the following lemma: 



Lemma 3.4. Let cr : R — > C &e a measurable function, continuous on R*, such that \a\ is locally 

2 ^ 
is locally integrahle at +00. Define i^l^) •= {'^ ^ ^'(R) : u = 



integrable at and that x 1— > 
T/; / e Li^o^(R) suc/i t/iat ||w 



1^ < +00} where <u,v>^ := LWiOr u{0 v iO d^- V C LUR), 



define ' " as the completion of for the norm \\ \\^. Then, the space (Lj.(R),<,>^ 



<.>o 



Hilbert space which also verifies £(R) ' " = L^(R). 

Proof. The fact that {L'^{M.), <,>^) is a Hilbert space is obvious. One needs only to show that the 
orthogonal space of f(R) for the norm || ||^ is equal to {0£(r)}. LetMinL^(R) be such that < u,v >^ — 

for all V in £(R). In particular, for all t in R, |cr(^)|^ u{^) il[o.t](0 d^ ~ 0. For all -0 in 



where if' denotes the derivative of V'. Thanks to the assumptions on \a\^ and x 1— > 
applies. Moreover, an integration by parts yields 



Fubini theorem 



m( f - e'^')dt)d^ = I |a(oi' m m d^ 



Thus < \a\^u,il} >= for all tp in ^(R). Since ^ ^ k('?)l^u(C) belongs to L[^^(R), it is easy to deduce 
that u is equal to 0. □ 



IS 



Remark 3.5. 1. Because the space ■y(R) is dense in L|^(R) for the norm \\ ||^ (see ISfp.lSj, it 
also possible to define the operator Mh on the space J'^(R) and extend it, by isometry, to all elements 
of L^{R). This is the approach of !ll5f and fS^ (with a different normalization constant). This clearly 
yields the same operator as the one defined by p.ip . However this approach does not lend itself to an 
extension to the case where the constant H is replaced by a function h, which is what we need for mBm. 
2. For the same reasons as in 1. it is possible to define the operator Mh on the space £{R) and extend 
it, by isometry, to all elements of (R) . Again, this extension coincide with p. II) . We will use this 
idea in sectionl 



11 



In view of p.3p . we find that 

< il[o,i],il[o,s] >if= — / :^:2H+i d^^RH[t,s). (3.11) 

As in the case of standard Brownian motion, one deduces that the process (_B'^) defined for all 

in M X 17 by: 

B^"Ht){u:) < c^, A/a(il[o,t]) >, (3.12) 

is a Gaussian process which admits, as the next computation shows, a continuous version noted _B'^^ := 
(i?'^) (i))jg^. Indeed, under the probability measure fi, the process is a fractional Brownian motion 
since we have, using p. lip and proposition 13.21 



E[B^"\t)B^"\s)]=E[< .,MH{il[o,t]) X .,A/h(i1[o,s]) >] 
=< Mij(il[o,t]), A/h(i1[o.s]) >l2(r) 

= < il[o,t], il[o,s] >H = Rnit, s). (3.13) 

Remark 3.6. The reason of the presence of the constant in formula p.ip is now clear since this 
constant ensures that, for all H in (0, 1), the process B'-"^ defined by (fX7g|) is a normalized fBm. 

Because our operator Mh is defined on a distribution space, we can not apply the considerations of |15j 
p.323ff about the links between the operator Mh and Riesz potential operator. However it is crucial for 
our purpose that Mh is bijective from L|f(]R) into L'^{M): 

Theorem 3.7 (properties of Mh)- 1- For all H in (0, 1), the operator Mh is bijective from L'j^(R.) 
into L^{R). 

2. For all H m (0,1) and {f,g) m {L^(R) n Ljj(R)f , 



< f,MH{g) >mm = < MH(/),.g >^,,^y (3.14) 



Moreover ( 3.14-\ remains true when f belongs to L]^^{S.) H L|^(M) and g belongs to (in 



this case ^3.14-\ reads < f^Mnig) >= < MH{f),g >l2(r); where <, > denotes the duality bracket 
between y"{R) and y{R) ). 

3. There exists a constant D such that, for every couple {H,k) in (0, 1) x N*, 
max|Mff(e,)(x)| < ^ (fc+1)'/^ 

Proof. 1. Since Mh is an isometry, we just have to establish the surjectivity of Mh, for all H in (0, 1). 
The case H = 1/2 being obvious, let us fix H in (0, l)\{l/2}, g in L^(R) and define the complex-valued 
function wfj on M by iu^(0 — -^=\£,\^^^^^ §{0 if C belongs to R* and u'|j(0) :~ 0. Define the tempered 

distribution vf^ by vf^ := ^Wh' where for all tempered distribution T, by definition, < T, f >:=< T, f > 
for all functions / in ^(K) and where f{x) = f{—x) for all x. We shall prove that vfj belongs to LfjCM.) 
and that Mh{v%) = .9- Note first that for all u in ^'(M), u = 27ru. It is clear that v^^ ~ w^^ and that 
v^jj belongs to Lj'Q^(M). Moreover, thanks to formula (|3.3p . we see that 

JR 

This shows that v^jj belongs to L|j(M). We can then compute Mh{vh) and obtain, for almost every ^ in 

R, 
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The previous equality shows that Mnivf^) is equal to g in and then establish the surjectivity of 

Mh. 

2. See equality (3.12) of [B]. The case where / belongs to Lj^^iR) n i|/(M) is obvious, in view of (P?Ti| . 
using the density of ^(R) in Ljj{R). 

3. is shown in lemma 4.1 of [15]. □ 

Of course if we just consider functions in L'jj{M.) instead of all elements of L'j^{R), the map Mh is not 
bijective any more. 

3.2 Study of ^ 



We now study the operator ^^jf- ■ It will prove instrumental in defining the integral with respect to mBm 
in section [5j 

Heuristically, we wish to differentiate with respect to H the expression in definition p.ip . i.e. differentiate 

the map H MH{u){y) on (0,1) for (w, y) in L|j(M) x M*, assuming this is possible. By doing so, we 
define a new operator, denoted from a certain subset of L'\j{R) to L^(M). Of course, in order to 

compute the derivative at ifo of iJ i—> MH{u){y), we need to consider a neighbourhood Vh,, of Ho in (0, 1) 
and thus consider only elements u which belong to P| L|^(M). However, as will become apparent, the 

formula giving the derivative makes sense without this restriction. 

In order to define in a rigorous manner the operator we shall proceed in a way analogous to the 

one that allowed to define Mh in the previous subsection. It will be shown in remark 13.91 that this 
construction effectively defines the derivative, in a certain sense, of the operator Mh. 
We will note c^j the derivative of the analytic map H ^ ch where ch has been defined in (|1.4p and set 

Ph := Let H belong to (0, 1). Define: 

- {u e ^'(M) -.u^Tf- / e LlM) and < +oo}, (3.15) 

where the norm || jU derives from the inner product on F/f (M) defined by 



1 



< u, V >s, := I Wh+ In \^\'^^" ^iO d^- (3-16) 

'^H Jr 

By slightly adapting lemma it is easy to check that (F/f (M), <, >^^(-]g-|) is a Hilbcrt space which 

verifies the equality Fjj(M) ^(M)^'^'-^ = £{M.)^'^'" . Note furthermore that, for every H in (0,1), 
the inclusion F^f (K) C L|^(R) holds. We may now define the operator ^^jf- from (F/f (M), <, >5„(r)) to 
(i^(R), <, >L2(jf)), in the Fourier domain, by: 



dM. 



OH 



f (u)(y) :== ~{Ph + In |y|) ^ \y\'^^'" u{y), for every y in R*. (3.17) 



In particular, one can check that, for / in .y(R), ^gjf iDiv) = ^-^ff (/)(2/) for almost every real y. 

Since ^1? (u) belongs to L^(R) for every u in F^f (R), is well defined and given by its inverse Fourier 
transform from (Fh(R), <, >5^(r)) to (^^(R), <, >i2(R)): 

aM^(u){x) = ^J-( ^^^'") )(-x), for almost every x in R. 



As in the previous subsection it will be useful to compute (/) for / in ^(R). We summarize, in 
following proposition, the main properties of 

Proposition 3.8. ^§jf- is an isometry from (THiR), <,>Sh{R)) i^'^i^)^ <^>l^{r)) which verifies: 



V/ e r«(R), = (3-18) 

y{f,g) e {rH{R)nL^R)f, < ^(/),g >^^^^^ = < /, ^(g) >^^^^^, (3.19) 

yf e S'{R)U£(R), and for a.e. xeR, ^if){x) ^ ^[MH{f){x)]. (3.20) 
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Proof. Equality p.lSp results immediately from the definition of ^^jf- and from (|3.17p . For any couple 
of functions {f,g) in (ThW H 



<^(/).3>z.^«) =2^<^(/).5> 



2ir^ an yyi ^L2(R) ^ ^ OH yyi ^l2(e)- 

It just remains to prove p.20p . Since we will not use p.20p in the sequel for / in £(M), we will just 
establish it here on ^(R). Let / be in ^(R) and H in (0, 1). Formulas ([321), (E21) and ([3^ can be 
summarized by 

MH{f){x) = 7_ff < \y\~'^'^^'^~^\ f{x + y) > for almost every real x, 

where we have written, by abuse of notation, \y\ '■^^^ for the tempered distribution y \y\ ^^^'^ 
and f{x + y) for the map y i— > f{x + y). Using theorem 15 . 1 71 which is given in subsection 15.51 below, we 
may write 

^[MH{f){x)]^iH < \y\~^'^^~"\f{^ + y) > +1H< \y\-^'/'~"hn\y\,f{x + y) > . (3.21) 
Furthermore, for almost every real x, we may write: 

= ^^¥ifn-^) = ^^y ^ -iPH + In \y\)^\y\^"'-"^I{y)){-x) 

= -PHMH{f){x) \^Hy ^ \y\^^'^'"^ In \y\Kym~x). 

Define, for every real x, I{—x) := F{y ^ \y\^^^'^~^^ \n\y\f{y)){—x) and, for every H in (0,1), vr '■= 
2 r'(3/2-i?)sin(f (l/2-i7)) + 7rr(3/2-i7)cos(f (l/2-i7). Using [5] (p.173-174) we get, after some 
computations, 

/(-x) = (-CH%/2^)7if < \yr^'^^'"^ In |?/|, /(x + ?;) > -^.^ < \yr^^^'-"\ fix + y) > . 
We finally have, for almost every real x. 



^-M^if)ix) = (Z2^ + ^2^) < \y\-^'/'-"\f{x + y) > < \y\-^'/'~"hn\y\Jix + y) > 



which is nothing but (13311) since ^22L£k ^ = 7'ff . □ 

Remark 3.9. For all positive real r and real H in (0, 1) let E^f ,.(M) be the space defined by 



T.H,r{^) := TH{R)f][ fl L|j,(M) p and I]h(R) := U J:H,r{'R) ■ It is possible to 

\H'elH-r,H+r] J re(0,min(H,l-H)) 

show that, for all f in Sfl-(R), (/)(■) ('resp. (/)(■) J ^■s ^'je derivative, in the L'^(R.)-sense, of 

Mnif) (resp. of Mnif) )■ Note moreover that the inclusion £(M) U ^(M) C f] ^h.t{^) holds. 

4 Wiener integral with respect to mBm on M 
4.1 Wiener integral with respect to fBm 

Similarly to what is performed in [TS] and [5] (in these works this is done only for functions of L|^(R)), 
it is now easy to define a Wiener integral with respect to fractional Brownian motion. Indeed, for any 



iNote that we have the equahty p| Lf, (M) = i^(K) n -Lg(M) for every [a, h] C (0, 1). 

H&[a,b] 
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element g in L|^(R), define {g) as the random variable < ., Mnig) >■ In other words, for all couples 
in Oxi|^(M): 

j''{g){Lj):^<u;,MH{g)>= I MnigK-s) dB{s)iu), (4.1) 

where the Brownian motion B has been defined just below formula (|2.5p . We call the random variable 
J^{g) the Wiener integral of g with respect to fBm. Once again, when is a tempered distribution 
which is not a function, g{s) does not have a meaning for a fixed real s and {g) is just a notation for 
the centered Gaussian random variable < .,MH{g) >■ 

Remark 4.1. Note that, for H in (0, 1), we are able to give a meaning to {g) only for elements g 
which belong to i|^(R) and not anymore for all elements g in J^'{M.) as was the case for J^g(s)dB[s) 
(see remark \2.4^ . 



4.2 Wiener integral with respect to mBm 

We now consider a fractional Brownian field A (A(i, H)^^^ i/)gMx(o i)' defined, for all {t, H) in Mx (0, 1) 

and all uj in il, by A{t,H){uj) := B^^\t^uj) < Mjj(il[o^t]) >. We wish to generalize the previous 
construction of the Wiener integral with respect to fBm to the case of mBm. This amounts to replacing 
the constant 7J by a continuous deterministic function h, ranging in (0, 1). More precisely, let Rh denote 
the covariance function of a normalized mBm with function h (see dcfinition ll.Sp . Define the bilinear form 
<, >f^ on X £(R) by < il[o,t] , il[o,s] >^ = Rh{t,s). Our construction of the integral of deterministic 
elements with respect to mBm requires that <, be an inner product: 

Proposition 4.2. <, is an inner product for every function h. 

Proof. Sec appendix [Bl □ 

Define the linear map Mh by: 

A'h: (£(M),<,>J ^ (L2(R),<,>^,(jj^) 

il[o,t] ^ A/,,(il[o,t]) := il4(t)(il[o,t]) := MHiil[oj])\H=h{t)- 

Define the process B^'^\t) = < ., Af/i(il[o^t]) >, i G M. As Kolmogorov's criterion and the proof of 
following lemma show, this process admits a continuous version which will be noted S'-'*^. A word on 
notation: we write S*- -* both for an fBm and an mBm. This should not cause any confusion since an 
fBm is just an mBm with constant h function. It will be clear from the context in the following whether 
the "ft," is constant or not. Note that: 

a.s., VteK, B('^)(t) = S(^)(t)|^^^^^^ (4.2) 
In view of point 2. in remark 13.51 we may state the following lemma. 
Lemma 4.3. (i) The process B^'*' is a normalized mBm. 

(ii) The map Mh is an isometry from (£(R),<,>^) to (L^(R), <, >^2(-r)). 
Proof. The process i?'-'*' is clearly a centered Gaussian process. Moreover, for all (s,t) in R^, 

E[B('')(t)S('')(s)] = E[< .,A4(il[o,t]) X .,Mhiil[o,s]) >] = < Ahiil^o^t]), Mh{il[o,s]) >^2(r) 

e**«)(l - e-^"«) 



— < M;,(t)(il[o,t]),M;,(,)(il[o.,]) >^,(^) - c,(,)c,(,) / TTa+^s 



Ch{t)Ch{s'\ 2 



By isometry, it is possible to extend Mh to the space £(R) and we shall still note Mh this extension. 
Define the isometry j'^ := C o Mh on f (R) ' '', i.e: 
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: (f(M)^'",<,>J (L^(R),<,>^.(R)) ^ ((L^),<,>(i.)) 

il[o,t] 1^ A4(il[o.t]) ^ < .,Mh(il[o,t]) > ■ 

Wc can now define the Wiener integral with respect to mBm in the natural following way: 

Definition 4.1. Let S^*) he a normalized multifractional Brownian motion. We call Wiener integral 
on M of an element u in £(R) with respect to B^'^^ the element j'^{u) of {Li^) defined thanks to the 
isometry j'^ given just above. 

Remark 4.4. It follows from definition |^. j| that the Wiener integral of a finite linear combination of 
functions il[o,t] is Q;fcil[o.tj.]) ~ "Y^^l^i akB[^\ Moreover, for any element u in £(R) 

(which may be a tempered distribution), the Wiener integral of u with respect to mBm, still denoted 

J^{u), is given by J^'^{u) ^= lim Sf'^{un), for any sequence of functions (u„) in £'(M)^ which 

n— ^+oo 

converges to u in the norm || ||^ and where the convergence of J^ {un) holds in (i^). 

Since we now have a construction of the Wiener integral with respect to mBm, it is natural to ask which 
functions admit such an integral. In particular, we do not know so far whether ^(M) C ^(M) . The 



next section contains more information about the space £{\ 



5 Stochastic integral with respect to mBm 
5.1 Fractional White Noise 

The following theorem will allow us to give a concrete example of a derivative of an (5)*-process. 

Theorem 5.1. 1. For any real H in {Q, I), the map Mh{iI[o,.]) : R ^ ^'(M) defined by MH{\l[o,.]){t) := 
M^f(il[o.t]) is differentiable ower R and its derivative, noted ^ [il/// (il[o,t] )] , is equal to 
s-k, where the convergence is in 

2. For any interval I of M. and any differentiable map F : I ^ J?^'(R), the element < .,F{t) > is a 
differentiable stochastic distribution process which satisfies the equality 
i<.,F{t) >=<.,-|F(i)>. 

Proof. The proof of point 1 is a mere re-writing of the one of lemma 2.15 of [4] by replacing by 
the operator Mh. Point 2 is theorem 2.17 of [3]. □ 

Let H e (0, 1). The process (B'^"'' (t))tm defined in (|XT^ is an fBm, and Af//(il[o,t]) belongs to L'^{R) 
for every real t. Hence, using equality (|3.14p . we may write, for every real t and almost surely: 



-t-oo 

B^"\t)^ <.,MH{il[o,t])> = < .,^< il/H(il[o,t]),efe >^,^g^ efc > 

fe=0 

+00 j-oo / \ 

= V< il[o.t],MH(efe) >^2/R> <-,efc> = V( / MHiek)iu)du) < .,ek> . (5.1) 
fe=o k=o V-^o / 

(|5.ip and the previous theorem lead to the definition of fractional white noise [111 15] : 

Example 5.2 (Fractional white noise ). Let: 

W^"\t) :=^ MnickM < .,ek >. (5.2) 

k=0 

Then (M^*-^-* (0)tgi{ a [S)* -process and is the (S)* -derivative of the process (-B*-^-* (t))^^^. 

The proof of this fact is simple: for any integer p ^ 2, using remark \2.3\. the mean value theorem and 

the dominated convergence theorem, 
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B(g)(t+,.)^_B(-»)(t) 



+00 



k=0 



t+r 



MH{ek){u) du~ MH{ek){t) < .,ek > 



J2 {- MH{ek){u)du-MH{ek)it)] 



+00 



-2p 



t+r 

{MH{ek){u) - M„{ekm) du ] ~> 0. 
r It I r-5.0 



(5.3) 



Remark 5.3. 

P ^ 2. 



In particular we see that for all {t, H) m M x (0, 1), W^'^'(t) belongs to {S-p) as soon as 



Remark 5.4. There are several constructions of fBm. In particular, operators different from Mh may 
he considered. uses an operator denoted on the grounds that fBm as defined here is not adapted 
to the filtration generated by the driving Brownian motion as soon as H 7^ 1/2. While this is indeed a 
drawback, the crucial property for our purpose is that the same probability space (J?^'(M), f/, /i) is used for 
all parameters H in (0, 1). This allows to consider simultaneously several fractional Brownian motions 
with H taking any value in (0, 1), which is necessary when one deals with mBm. We choose here to work 
with Mh rather than with and of 2] as its use is simpler. and would nevertheless 

allow for a more general approach encompassing the whole family of mBm at once. This topic will be 
treated in a forthcoming paper. 



5.2 Multifractional White Noise 

The main idea for defining a stochastic integral with respect to mBm is similar to the one used for fBm. 
We will relate the process B^''^ to Brownian motion via the family of operators {Mh)^u£{o i))- This 
will allow to define a multifractional white noise, analogous to the fractional white noise of examDle l5.2l 
From a heuristic point of view, multifractional white noise is obtained by differentiating with respect to 
t the fractional Brownian field A{t,II) (defined at the beginning of section |4?2|) along a curve {t, h{t)). 
Assuming that we may differentiate in the sense of (5)* (this will be justified below), the differential of 
A reads: 



dA{t, H) = |^(t, H) dt + ^{t, H) dH = -^{t) dt + |^(t, H) dH 

= W^"\t)dt+^{t,H)dH, (5.4) 

where the equality will be shown to hold in (S)* . With a differentiable function h in place of H, this 
formally yields 

dA{t,h{t)) = (W('^«)(0 + /i'(t) ^(i,i?)|«=,(,)) dt. (5.5) 

In view of the definition of the stochastic integral with respect to fBm given in [15] , [H] and [5] , it then 
seems natural to set the following definition for the stochastic integral with respect to mBm of a Hida 
process X : K (5)*: 



[ A(s)dB('')(s) [ X{s)dA{s,h{s)) 
Jr Jr 



X{s) o (Ty(''(^))(s) + h'{s)^{s,H)\H^,^,~,)ds. (5.6) 

We shall then say that the process X is integrable with respect to mBm if the right hand side of (|5.6p 
exists in (S)* . Remark that when the function h is constant we recover of course the integral with respect 
to fBm. 
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In order to make the above ideas rigorous, we start by writing the chaos expansion of B^'^\ Since Mnig) 
belongs to L^{R) for all (.g, H) in y{R) x (0, 1), we may define, for aU H in (0, 1), Mh : ^'(K) J?^'(]R), 

by 

<MH(uj),g> = <u,MH{g)>, for ^ - a.e. w in f7 = .J^'(M). (5.7) 

Moreover, in view of remark [5751 we may extend (|5.7p to the case where g belongs to L^^ (R) by writing, 
for all g in L|^(E) and almost every w in £7, 

< MH{io),g>:^ lim < Mi/(a;),.g„ > = lim <uj,MH{gn)> = <uj,MH{g)>, (5.8) 

for every sequence ((7n)„gpj of functions of oS^(M) which converges to g in the norm || ||^2 (g-,. For all real 
t and integer k in N, define the element of L^^s^{M.): 

4^ ■■= ^i^i^k)- (5-9) 

It is clear that, for all t in R, the family of functions [df^)^^^^ forms an orthonormal basis of L^jj-|(M). 
Let us now write the chaos decomposition of mBm. For almost every uj and every real t we get, using 
(Pl^ and (jO)) . 



> 



+ 00 



4-00 



+ CXJ 

k=0 

+ 00 

= ^< A4(t)(il[o,t]),efe >^ < w,A4(t)(4*^) > . 



fc=0 

We get finally: 



+ 00 / „t 

.s, ViGR, i?(''Hi) = E / ^^^'^ 



(t)(efc)(s)c?s < .,efe > 



(5.10) 



We would then like to define multifractional white noise as the (5* )-derivative of B'-'^\ which would be 
formally defined by: 



+ 00 



fe=0 



(5.11) 



assuming h is differentiable. The following theorem states that, for all real t, the right hand side of (|5.1ip 
does indeed belong to (S)* and is exactly the (5)* -derivative of B^^^^ at t. 

Theorem-Definition 5.1. Let ft, : R -> (0, 1) be a function such that the derivative function h' is 
bounded. The process lyCO := (W^('')(0)t6R defined by formula (TOTIl is an (S)* -process which verifies 
the following equality in (S)* : 



+00 +00 , „t 

M^W(i) = E^^M*)(efc)W<->efc>+ft'(i) E / 
Moreover the process B'--^^ is {S)* -differentiable on R and verifies in {S)* 

^(" = »'"-'<" = sW'' "<'»!■ 

In order to prove this theorem, we will need two lemmas. 



< .,6^ >. 



(5.12) 



(5.13) 
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Lemma 5.5. For H m (0, 1) and f in ^{R), define 5/ : M x (0, 1) K by gf{t, H) := [ MH{f){x)dx. 

JO 

Then 

(i) The function gf belongs toC°°(Rx (0, 1), R), 

r+oo 

(ii) Vx e M, Mh{I){x) ^an u"-^l'^{f{x + u) - fix - u)) du. 

Jo 

where an has been defined by ^3.10^ . 

In particular, the function {x,H) Mi{{f){x) is differentiable on R x (0, 1). 
(Hi) Assume that : R — > (0, 1) is differentiable. Then, for any real tg 



d 



to 



j^[gf{tMml^,,=Muito){f){to) + h'{t^)j ^{f){s)\^^^^^^^ds. (5.14) 

Proof. (i) Define ^/ on R x (0,1) by fJ.f{t,H) := J^{u - t)\u - tf^^^^ f{u) du, for / in ^(R). 
Using ((X^ we get, for all (t, il) in R x (0, 1), the equality 

gf{t, H) = a^^[/i/(t, H) - ;U/(0, H)]. (5.15) 



A change of variables yields 



/•+00 

^if{x, H) = - I \t~ xf-'^'f{t)dt + I \t- x\"-'/'f{t) dt 

J —00 

,H-l/2. 



u^-^l\f{x^u)- f(x-u))du. (5.16) 







Thanks to (|5.15p and to the fact that the map ?/ 1— >■ Ofj, is C°° on (0, 1), it is sufficient to show that 
the function /i/ belongs to C°°(R x (0, 1)). In view of applying the theorem of differentiation under the 
integral sign, define j{x, H,u) := u^~^/^{f{x + u) — f{x~u)) for u in R^. Let n in N and (ai, 0^2) in 
such that ai + a2 = n. For almost every u in R!i_, {x, H) 1-^ j{x, H, u) is C" on R x (0, 1) with partial 
derivatives given by 

{X,H,U) = (lnw)"^U^^l/2/j(ai)/ ^ X j(ai)( 

Fix {xq, Hq) in R X (0, 1). Let us show that /i/ is C" in a neighbourhood of {xo,Hq). Choose (a, b) such 
that a < Xq < b and Hi,H2 such that < Hi < H < H2 < 1. We have 



I ^ u „M ^ u,i^^i-i/2 

' dx"^dH°"^ 



-{x,H,u)\^\uf'-'^^\\nur,l^o<u<i} sup \f('^^\x±u)\ 

ix,u)e[a,b]x[0,l] 



+ |«||ln?.r^|/("i)(x±u)|il{i^l„|}, (5.17) 



where /'^"^^(x ± u) := /("^^(.t + u) ~ f^°''^\x — u). A Taylor expansion shows that there exists a real 
D such that, for all {u,x) in R x (a, 6), \u\ \f^°'^\u±x)\ ^ D. As a consequence, there exists a real 
constant C such that, for almost every u in R^ and every {x, H) in [a, b] x [i/i, iJ2], 



3"j 



(x,g,^)| sc: C [l^^l^^-i/^llni^r^ il{o<„<i} + llni^r^ ^ il{i^u}]. (5.18) 

Since the right hand side of the previous inequality belongs to L^(R), the theorem of differentiation 
under the integral sign can be applied to conclude that the function /j,/ is of class C" in [a, b] x [Hi,H2], 
for all integer n and all / in This entails (i). 

{a) and yield 

d f~^°° 
MH{f){x) = aH^[tif{x,H)]^aH j u"-^'^{f' {x + u) - f {x - u))du, (5.19) 

which establishes (ii) and the fact that {x,H) MH{f){x) belongs to C°°(R x (0, 1)). 



19 



(Hi) For a differentiable function h, we have, for every real to, 

|[5/(^,Mi))]L.,„ = ^(^o,Mio)) + /i'(to)||(to,Mio)). 
((5A5)) . ((5l9)) and (|3:20)) show that, for every / in ^(M), 



and finally: 



d 







to 



□ 



Lemma 5.6. The following inequalities hold: 

(t) y[a,b] C (0,1), 3peM.: V/c G N, sup ^{ek){u)i^p{k + lf^^ln{k + l). 



(ii) yt e K, Vr e M;, BAW G M, Vfc G N : sup \£[ge^{u,h{u))]\ AW (fc + if^^- 

ue[t-r,t+r] 

Proof, (i) Since is a subset of ThC^), (|XT71) entails that ^^(e^) belongs to n L2(E) for 

every k in N. Furthermore ek{y) = {—i)'^ ^ \/^e-k{y) for every integer k in N* and for almost every real 
y (see lemma 1.1.3 p. 5 of [37])- Thus, for every H E [a, b] and almost every u G R, 

^^(«)=2^^^^(-^) = -27 / e'-HPH + ln\y\)^\y\'^'""€k{y)dy 



^ / e^-yMH{ek)iy)dy ~ ^ / e^-y\y\''^-"{\n\y\){^^f-'ek{y)dy 

-(3h MH{ek)iu) + / e^'^y\y\'/'-"{\n\y\) e^iy) dy . (5.20) 



= .Vk{u) 

Then, using (j2.8p . we see that there exists a family of real constants denoted (/50i<i<ii such that we 
have, for every couple {H, k) in [a, &] x N and almost every real w, 



^ Pi [S\y\^2Vk+i \y\"^'"\^^ \y\\\eu{y)\dy + 4|>2yFFT \y\"^~" i^^ \y\)\eu{y)\dy 

r (.2Vfc+l p+oo 



(5.21) 



An integration by parts yields 

=p3((p4+lnfc)(l + fc3/4-^/2))^p5fc3/4-ff/2lnfc. (5.22) 

Using the change of variables u = 2/ a/T 7 we get /j''^ ^ Pe |wP^^~^(lnu)e~"^dM = pe "^2^/!"^' where 

j]"'' := /^^°° |w|"lnti e^" du. When 5 > 3e, an integration by parts shows that < ^""^6^"^ InS 
for < a < 1/2, and that j]"' ^ 5" e""' Inu du for -1/2 < a < 0. Hence we get 

^ |^i/4_^/2 for l/2<i7<l ^ ^ 
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and we finally obtain 



In k 



^2 ^Ps (5-24) 

Using (j5.20p to (j5.24p . item 3 of theorem l3.7l and the fact that both functions H i-^ (3h and H ^ are 
continuous on [a, h] we get, for every k in N, 



sup 

(i/,ii)G[a,6]xI[ 



dM. 



-mriekKu) ^P9ik + if^ + pio sup [((fc + i)'/^-^/V(fc + 1)1/^-^/2) hi(fc + 1)1 

He[a,b] ^ 

(fc + l)^/^ln(fc + l). (5.25) 
(m) Let i e M and ?■ > be fixed and define [a,b] :— [ inf h{u), sup h{u)]. Using (|5.14p we 

ue[t-r,t+,-] u£[t-r,t+r] 

have, for every k in N, 



sup \^[9e^{u,h{u))]\ ^ sup \MH{ek){u)\ 

u£[t-r,t+r] (H,u)£[a,b]x'R 

+ i\t\+r) sup \h'{u)\ sup 

uelt~r,t+r] (ff,ti)6[a,b]xl 

The result then follows from (i) above and item 3 of theorem 13.71 □ 

2/3 

Remark 5.7. In the sequel, we will only need the bounds of (i) and (ii) with {k + 1) in lieu of (fc + 1) . 
We may now proceed to the proof of theorem 15. II 

Proof of theorem 15.11 From equality (|5.1ip defining multifractional white noise and equality (|5.14p . 
we have formally 

+00 H-oo „i 

W^('')(i)=5^M,(,)(efc)(i)<-,e..> + /»'WE / ^T^(^)lif=Mt)^^<-'^'^>- (^.26) 

fc=0 fc=0 ° 

In order to establish that W^'^''{t) is well defined in (S)* and that equality (|5.12p holds, it is sufficient to 
show that both members on the right hand side of the previous equality are in (5)*. 

-|-oo 

For t in M, definition (j5.2p of fractional white noise shows that ^ M,j(t)(efc)(i) < .,6^ >= M^^'''*^^(i) 
and thus belongs to (5)*. 



k=0 



Let us show that Vff(t) := ^ / ^^^^^^{s)ds < .,ek > belongs to (i^). Using (j3.19l) we may write 



fc=0 ' 



ES<e.,^(il[o,])> 



LH 



IIw"(i1m)III2(r) - l|il[o.t]ll5^ 



< +0O. 



As a consequence. W'^^^t) is the sum of an (5)* process and an (L^) process, and thus belongs to {S)* . 
Wc are left with proving equality ((gl^ . i.e. that W'^^^t) is indeed the (5)* derivative of S('')(t) for 
any real t. 

Let r ^ and < ^ (the case t < Q follows in a similar way). The equality l^('')(i) H^(''(*))(t) + 
/i'(i) V/i(t)(t) and remark [5731 entail that VF'^''-'(t) belongs to [S-p) as soon as p > 2. For such a p, one 
computes: 



5e,(t + r,/i(t + r)) -<7e,(i,/i(t)) 



+00 

E 

+ CX3 



< .,ek > 



1 



2p 



- |K(t,Mt))] 



(5.27) 
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Using lemma [^751 and the Mean- Value theorem we obtain, for r in (—1/2, l/2)\{0}: 

Jp,rm ^(SFTSj^fs sup \£[g^^{uMu))]\\ ^'^D^^ 

^ ' \ Mg[t-l/2,t+l/2l / \ ^ J 



22(P-1) (fe + l)2{p-l) ' 

where D := Dt{l/2). Since Jp^r k{t) — > 0, equality (|5.13p follows from the dominated convergence 
theorem. 

Remark 5.8. In (ii) of lemma [576[ the real constant Dt{r) can be taken independent oft if the function 
t t-^ t h'{t) is bounded over M. 

We note that multifractional white noise is a sum of two terms: a fractional white noise that belong to 
(5_p) as soon as p ^ 2, and a "smooth" term which corresponds to the derivative in the "77" direction. 
This is a direct consequence of the fact that the fractional Brownian field A(t, H) is not differentiable in 
the t direction (in the classical sense) but infinitely smooth in the H direction. 

Proposition 5.9. For p ^ 2, the map 1 1-> ||W^^''HOII-p continuous. 

Proof. By definition. 



Using the estimate given in lemma 1^751 (ii), we sec that \\W'^''Ht)\\_p is the sum of a scries of continuous 
functions that converges normally on any compact. □ 

5.3 Generalized functionals of mBm 

In the next section, we will derive various Ito formulas for the integral with respect to mBm. It will be 
useful to obtain such formula for tempered distributions. In that view, we define generalized functionals 
of mBm as in 0]. 

Theorem-Definition 5.2. Let F be a tempered distribution. Fort in R!j_, define 

FiB'^'^Ht)) := -j=— ^ (fc!)-it-2^"(*) < F,C(*.,(*).,) > h (( MMt)(il[o,*]) f") (5.28) 



where the functions £,t,H,k o,re defined for {x,H,k) m M x (0, 1) x N by 

6,ff.fc(x) := iV2f' t'^" hu{x/{V2t")) cxp{-^} = 7:''\k\fh^"cM-^}eu{x/{V2t")). 

(5.29) 

Then for all t inWj_, F{B^[t)) is a Hida distribution, called generalized functional of B^^\t). 

Proof. This is an immediate consequence of [23 p. 61-64 by taking / := i\//i(il[o.t]). □ 

Remark 5.10. As shown in when F = f is of function type, F{B^{t)) coincides with f{B''^{t)). 

The following theorem yields an estimate of \ \F{B^''''' {t)\f_^ which will be useful in the sequel. 

Theorem 5.11. Let h : R ^ [Hi,H2] C (0,1) be a continuous function, B'--^^ an mBm, p G N and 
F S J?fp(M). Then there is a constant Cp ^' ^ , independent of F, such that 

yt > 0, \\FiB^^Ht))\\'_j^ s; max{t-4p''W,i4pMt)}i-Mt) cj^HuH,)^p^2^^^ (5_3q) 
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Proof. For H G (0, 1) and p € N, Theorem 3.3 p. 92 of [3] ensures that there exists Cp such that, 
Vt > 0, 

Now if H belongs to \Hi,H2], it is easy to show, by examining closely the iteration of (23) p. 94 in [3], 
that one can choose a constant C^^^'^^^ independent of H. We hence have 

Vt>0, yHe[Hi,H2], \\F(B"{t))\tpi^max{t-^P",t^P''"}t-" C^^^'^^^-^lFtp. (5.31) 

For t > 0, one only needs to set H = h{t) in ((OT|) to get ((00)) . □ 

5.4 S-Transform of mBm and multifractional white noise 

The following theorem makes explicit the 5'-transforms of niBni, multifractional white noise and gener- 
alized functionals of mBm. 

We denote by 7 the heat kernel density on IR+ x R i.e ■f{t, x) := ^^^^ cxp {^^^} if t 7^ and if < = 0. 

Theorem 5.12. Leth:M.^ (0, 1) he a function and (SC') (i))jgjj (resp. {W'^^\t))^^^) he an mBm 
(resp. multifractional white noise). For rj G ^(K) and t £ M, 

(i) S[B''''-^t)]{i]) = < 1], AIh{il[o,t]) >£,2(j{') — gri{t,h{t)), where g,^ has heen defined in lemma [57Si 

(u) 5[W^('')(t)](r?) = f^[g^{t,hm - Ah^tMit) + h'{t)J^ 9M^{r^)(s)\^^,^^^^ds. 

(Hi) ForpeN and F e ,yLp(M), 

SiFiBC^HtMTj) = (f,^ (^t^'^('\.- A4(,)(,y) (u) du 

Furthermore, there exists a constant Cp, independent of F,t and rj, such that 

ISlFiB^f'^tmrj)]^ ^ max{t-4p''(*),t4ph(t)|^-h(t) j^j^^ cxp{\APrj\l}. (5.32) 

Proof, (z) Thanks to (|5.10p and lemma [^31 we have, for every rj in ^(M) and t in R, 
+00 

S{B^''Ht))iv) ^4(t)(il[o,t]),efc >i2(K)< V,ek >mR) = < A4(t)(il[o,t]), ry >^2(k) 

fe=o 

= < ^'^[o.t],Mm)iv) >i2(R) = .9r,(i,^(0)- 
(m) This is a straightforward consequence of lemma (j5.13p and (i). 

(Hi) The first equality results from theorem 7.3 p. 63 in 23 with / = Af^(t)(il[o,t]) and from the equality 
(i). Equality ([5732]) rcsuhs from ((O0| as in theorem 3.8 p. 95 of [4]. □ 

Remark 5.13. Using lemma W^ and \5.11^ we may also write: 

V(t,77)eRx^(M), 5(W^W(t))(,,)^^ ^K(t,Mi))] <^?,efc>^2(R). (5.33) 

fc=0 

5.5 The multifractional Wick-Ito integral 

Wc arc now able to define the Multifractional Wick-Ito integral, in a way analogous to the definition of 
the fractional Wick-Ito integral. In the sequel of this work, we will always assume that h is a function 
on R with hounded derivative. 
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Definition 5.1 (The multifractional Wick-Ito integral). Let y : M — > (S)* be a process such that the 
process t i-> Y(t) o W^'^^t) is (S)* -integrable on R. We then say that the process Y is dB^'^^ -integrable 
on K or integrable on M with respect to mBm i?^'') . The integral of Y with respect to B^'^'^ is defined by 

Y{s) dB^^\s) := [ Y{s)oW^''^s)ds. (5.34) 
For a Borel set I ofR, define j ^Y {s)dB^^\s) il/(s) Y {s)dB'^^\s) . 

When the function h is constant, the multifractional Wick-Ito integral coincides with the fractional Ito 
integral defined in [TS], [5], H] and [5]. In particular, when the function h is identically 1/2, (|5.34p 
is nothing but the classical Ito integral with respect to Brownian motion, provided of course Y is Ito- 
integrable. The multifractional Wick-Ito integral verifies the following properties: 

Proposition 5.14. (i) Let {a,b) in a<b. Then jj^ dB'^^\u) = B^'^^b) - B^''\a) almost surely. 

(a) Let X : L ^ (S*) be a dB^'^^ -integrable process over L , a Borel subset of W. Assume Jj- X(s)dB'^'^\s) 
belongs to (L^). Then E[/^ X(s)dB('') (s)] = 0. 

Proof, (i) From (m) of theoremEU 1 1-^ S{il[a,b]it) W^''^ {t)){T]) is measurable on R for any r] in ^(M). 
Moreover, for any integer po ^ 2, we have 

|5(il[,,,](t) W^'^Hmv)] ^ l|M^('Ht)ll-poe*l< 

thanks to lemma [275] By proposition 15.91 t !-> llW^'-'^HOIUpo continuous thus integrable on [a,b]. 

Theorem 12.91 then entails that t H> il[a,6](i) W^'^^t) is (iS*)-integrable over M. It is easily seen that the 

S-transforms of dB^^'^^u) and B^^^{b) — B'^'^^a) coincide. The result then follows from the injectivity 
of the S-transform. 

(a) One computes: =^[go(s,^(s))]=o 

S{f X(s)dB('')(s))(0) = / S{X{s)){0) 's{W^''Hs)){Qj ds = 0. 



Now, when Jj X{s)dB^''^s) belongs to (L^), E[J^ X {s)dB^''^s)] = S{Jj X{s)dB^'''>{s)){0). □ 

Theorem 5.15. Let L be a compact subset o/M and X : 1 1-^ X{t) be a process from I to (L^) such that 
t I— > S{X{t)){'if} is measurable on I for all rj in S^{M.) and t i— > ||X(t)||g belongs to L^{I). Then X is 
dB^'^") -integrable on I and there exist a natural integer q and a constant Cj such that, 



X{t) dB'^''\t) 



i^Ci / \\X{t)\\^dt. (5.35) 



Proof. For 77 e ,-^(M), the measurability on / of t S{X{t) o W^^\t)){vi) is clear since 
S{X{t) o W^^\t)){ri) = S{X{t)){ri) j-^[g,^{t, h{t))]. By lemma we have, for any integer g ^ 2, 

\s{x{t) o w^'^Hmv)] < ii^wiio m^'^Hm-, ei"i? 

for every t in /. Since t H' ||M^^''^(i)||_q is continuous by proposition 15.91 and t i-> ||X(^)||p belongs to 
L^{I) by assumption, the result follows from theorem l2.9l We refer to theorem 13.5 of [57] for the upper 
bound. □ 

Remark 5.16. One can show, using avvendix m that inequality l\5.35^ is true for every integer q ^ 2. 

It is of interest to have also a criterion of integrability for generalized functionals of mBm. In that view, 
we set up the following notation: for p g N, < a < &, we consider a map F : [a, b] — > ^_p(R) (hence 
F{t) is a tempered distribution for all t). We then define F{t,B'^^\t)) := F{t){B^^\t)). RecaU the 
following theorem (see [T0|, lemma 1 and 2 p. 73-74): 
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Theorem 5.17. Let I be an interval of M., t i-> F{t) be a map from I into ^_p(R), t H> (p{t,.) be a 
map from I into .y{R) and tg € I. If both maps t > F{t) and t i— >■ ip{t, .) are continuous (respectively 
differentiable) at to, then the function t i— >< F{t),(p(t, .) > is continuous (respectively differentiable) at 

to. 

Theorem 5.18. Let pSN, 0<a<6 and let F : [a,b] — )• S^^p{M.) be a continuous map. Then the 
stochastic distribution process F{t, B'^'^'^ (t)) is both (S)* — integrable and dB^^^ -integrable over [a,b]. 



Proof. Wc shall apply theorem 12.91 

The measurability of 1 1-> S[F{t, B^'^^ {t))]{ri) results from {Hi) of theorem 15. 121 the continuity of the two 
maps t K-s- F{t) and t i-s- 7 (t'^'^^^\ . — J* Mm^{7]) (u) du^ and theorem [5. 171 
Since h is bounded on [a, b], lemma [231 and (fOO)) yield 

s;max{t-2p''(*),<2M*)}i-M*)/2 ^/^f^ max|F(s)|_^ exp{i|^P77|^} 

s^[a,b] ^ Z 

s^[a,b\ Z 

(5.36) 

where Hi := min h(s) and H2 '.= max /i(s). This yields the second condition of theorem l2.9l and shows 

s^[a,b] s^[a^b] 

that F{t,B^^\t)) is (5)*-integrable over [a,b]. 

For (i_B'^'')-integrability, we first note that, by theorem 15. 121 (ii). 

S[F{t,B^'^Ht))oW('^Ht)M ^ S[F{t,B^'^Htmv) f^bnit^m)]- 

Since the function i ^ [9r]{t, h{t))] is continuous (by lemma 1575)) . the measurability of 1 1—> S[F{t, B'-'^\t))o 
VF^''-'(t)](77) for every function 77 in J^{M.) follows. 

Moreover, for every integer po ^ max{p,2}, F{t, B'-'^'> (t)) and W'-'^')(t) belong to (iS_p„) for all t in [a,b]. 
Using lemma [^751 and (j5.30p . we may write, for all t in [a, 6], 

\s[F{t,B^>^Ht))ow^'^Htm\ ^mt,Bwm\^^^ m^'^Hm-p, e^piivO 



where K sup HW^e^HOILp^ {{^f'""' + b''""^m)^ + b"^) ^^P^ max |F(5)|_^^. 

te[a,b] se[a,b] 

Theorem 12.91 applies again and shows that t F{t, B^'^\t)) o W^^\t) is integrable over [a, b]. □ 

Remark 5.19. Recall that a function / : M — > K is said to be of polynomial growth if there is an integer 
m inN and a constant C such that for all x G M, |/(a;)| ^ C(l + The previous theorem entails in 

particular that both quantities f{B'^''\t)) dt and £ f{B^^'\t)) dB^^'\t) exist in (S)* if f is a function 
of polynomial growth. 



Example 5.20 (Computation of B^''^(t) dB^^^t)). Let T > fixed. Then 

(5.38) 



B^''\t)dB^''\t)^ r W^'^Ht) o B'^''\t) dt ^ r ^^^^ o B^''\t) dt. 



dt 

Let us prove that the last quantity is equal to ^B'^'^^Tf^ := i {B^''^T)oB^''^T)) ^ \ (SC^' (r)^-r2/i(T)-) 
(see remark \2.'l\ . It is sufficient to compute the S-transforms of both sides of the equality. For rj in o5^(R), 

^ B^'^\t)dB('^\t)j (r^) ^ ^ SiB^'^Hmri)SiW<^''Hmri)dt = ^ g,{t,h{t))i[g,it^hmdt 
= \gHTMT)) = i(5(i?(")(T))(,y))' = \s(b('^\T)oB^^\T)) (rj) 

= 5 Q (i3W(r)' - r2h(T))^ (^^^y (5 39) 
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To end this section, we present a simple but classical stochastic differential equation in the frame of 
mBm. 

Example 5.21 (The multifractional Wick exponential). Following \15\ formula (4.8) and example 
3.6, let us consider the multifractional stochastic differential equation 



{ dX{t) - a{t)X{t)dt + mrndB^"^^ it) 
I ^(0) e (S*), 

where t belongs to M_|_ and where a : M — > R and /? : M — > M are two deterministic continuous functions. 
( |5.^0P is a shorthand notation for 

Xit) = X{0)+ f a{s)X{s)ds+f /3(s) X(s) dB('')(s), (5.41) 







where the previous equality holds in (S) . Rewrite the previous equation in terms of derivatives in (S) 



as: 



^(t) = ait) X{t) + m X{t) o it) = iait) + mW^'^^ it)) o X(t) 

XiO) e is*)- 

We thus are looking for an iS*) -process, noted Z, defined on such that Z is differentiable on 
and verifies equation \5.42\ in (5*). As in it is easy to guess the solution of \5.4'/^ if we replace 
Wick products o by ordinary products. Once we have a solution of ( |5.^i^p , we replace ordinary products 
by Wick products. This heuristic reasoning leads to defining the process Z := (^(i)tgR^) by 

Z(t) :=X(0)oexp* (^^ a(s)ds + ^ ^(s)dB('')(s)^, t € R+, (5.43) 
where exp* has been defined in section \2.4\ 

Theorem 5.22. The process Z defined by \5.Ji-S^ is the unique solution in (5*) of l\5.42\ . 

Proof. This is a straightforward application of theorem 3.1.2 in [24]. □ 



Remark 5.23. (i) \24i uses the Hermite transform in order to establish the theorem. However it is 
possible to start from {5.41^ , take S -transforms of both sides and solve the resulting ordinary stochastic 
differential equation. 

iii) Equation \5.4(][ may be solved with other assumptions on a, /3. We refer to \29\ for more on stochastic 
differential equations driven by mBm. 

Remark 5.24. In particular when XiO) is deterministic, equal to x, aQ = a and /3() = /3 are constant 
functions, the solution X of l\5.41i reads 

X(t)=xexp{/35(^)(t)+ai-i/32t2'«(t)}^ t € R+, (5.44) 
which is analogous to formula (3.31) given in [Q in the case of the fractional Brownian motion. 

5.6 Multifractional Wick-Ito integral of deterministic elements versus Wiener 
integral with respect to mBm 

In section l¥^ we have defined a Wiener integral with respect to niBni. It is natural to check whether this 
definition is consistent with the multifractional Wick-Ito integral when the integrand is deterministic. 
More precisely, we wish to verify that J^fis) o W is) ds = J^if) for all functions / such that both 
members of the previous equality exist and that the left-hand side member is in (L^). In that view we 
first prove the following theorem. 
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Theorem 5.25. Let / : M — > R 6e a deterministic function which belongs to Lj^^(R). Let Z :~ {Z{t)^^-^ 

be the process defined on M by Z{t) := J{s)dB^''\s). Then Z is an {S*) -process which verifies the 
following equality in {S*) 



/ f{s)dB('^Hs)=Y^( . 



fis)f^[g,,{s,his))]dsj<.,ek>. (5.45) 

Moreover Z is a (centered) Gaussian process if and only if '^t^ (/g /(*) Jjlffefc (s, ^(s))]'^s) < +oo, 
for all t. In this case, 



Z(t) = f{s)dB^'^\s) ^ M WY.\J^ f{s)£[geds,h{s))]dsj 1 , Vi G M. (5.46) 

In particular, the process Z is Gaussian when the function f belongs to (R, R) and is such that 
sup|/'(0| < 

Proof. We treat only the case t e R^. The other case follows shnilarly. Let / be in Lj'Q^(R). In order 
to show (|5.45p let us establish a), b) and c) below, 
a) s H> /(s) o W'^'^\s) is (5)*-integrable over [0,t]. 
Let 77 € ^(R) and s in [0,t], using lemma 1^31 we get: 



\S{fis)oW^'^Hs)m^\f{s)\ SiW^'^H-^v) ^|/(.s)| \\W^''Hm-r,^e' 



2 I'd 



for s in [0,t] and for pq ^ 2. Since L is the product of a continuous function and a function of _Lj'Q^(R), 

a) is a consequence of theorem 12.91 

b) */ E£o (/o /(s)s[ffefc(s,/i(s))]'^s) < -,6*; > belongs to (5_pJ as soon as po ^ 2. 
Lemma inSl entails that there exists a real D such that, for every po ^ 2, we have 



Jo ■^(") d^jg-J^.M^))] rf^ \ ^||.||2 .g[o,ti 

(2fe+2)P0 I =^ II./ llLi([0,t]) Z^fc=0 (2fe+2)^fo 

c) /J/(.)dS('')(.) = ES (/^*/(,)^[5^J,,/,(,))]ds) < > in (S*). 

Denote := /J /(s)dB('') (s) = (ES /(s)^ [Se. (s, /^(s))]< -,6^ >) ds and define the (5*)-process 
T : [0,t] ^ (5*) by t{s) := /(s)ilK (s, /i(s))] < -.efe >. Moreover, for N in N*, define on [0,t], 

TN : s ^ tn{s) := J2k=o fi^)jE^9ek{s,h{s))]< .,efc >. Obviously we have, in (5*), $/ = T(s)ds, 
= lim J* T'jv(s) (is. It then remains to show that 'I'/ = lim J* T]y{s) ds in (5*). Let us use, for 

this purpose, theorem lA. 21 Let po be an integer greater than or equal to 2. It is easily seen that t„ and 
T are weakly measurable on [0,t] for every n in N (see definition lA.ip and that, t„(s) and t(s) belongs 
to {S-p„) for every n in N and s in [0,t]. Moreover, both functions s ||7Vi(s)||_pj^ and s i-> ||t(s)||_p^ 

belong to L^{[0,t],du) since \K{s)\\_^^ Ms)\\-p, ^ 1/(^)1 VeS (2^ + 2)-2(po-i) for a certain D 
given by lemma [^751 fw). We hence have shown that both functions t„(.) and r(.) are Bochner intcgrable 
on [0, t]. Besides, for every (n, m) in with n ^ m, we have 
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|t„(s) - r„i(s) 



ds = 



J2 /(«)il[5ejs,/i(s))]< .,efe > 

k—ra-\-l 



ds 



-Po 



+ CXD 



E /(s)i^K(s,Ms))]<.,e,> 

k—m-\-l 



ds < M 



+00 sup 1^ 



|[5e,(s,/l(s) 



1/2 



(2fc + 2) 



2po 



+ CXD 

E 



1 



\ 1/2 



1 {2k + 2) 



2(po-l) 



{ ?i , m )—)■(+ cxj ,+ 00 ) 



^0, 



where il/ II/IIli([o t]) ^^'^ ^ again given by (ii) of lemma [5.61 Theorem IA.2I then applies and 
allows to write that lim f* tm(s) ds = f* t(s) ds in (S*). We hence have shown that ^/ = 

lim tn{s) ds ~ jQT{s)ds = in {S*). This ends the proof of c) and establishes formula (j5.45p . 

2 

< +00, for all t, then Z{t) is the (L^)-limit of a sequence of 
independent Gaussian variables. Formula (|5.46p is then obvious. When / is of class and such that 

+00 / «y- \ 2 



If Etro{lofi^)Tsi9e.{s,his))]ds 



sup|/'WI < 



+00 / „t 

-00, an integration by parts yields that I / /(s)^L9' 

fc=o V-'o 



Xs,h{s))]di 



< +00. 



□ 



It is easy to check that definitions 14.11 and 15 . 1 1 coincide on the space £(R). Indeed for / := 5I]fc=i '^k'^^lo.tk] 
in £(M), remark H31 and equality (|5.13p entail that J^^{f) = 'Y%=iOLkB\^^ almost surely. According to 
(i) of proposition 15.141 we have the equality f {s)dB^''\s)ds = Y^^^-^^ak Jjg,il[Q^tk]is)W^'^'> {s)ds = 
X]fc=i "^fc^tl'^ almost surely. This implies in particular that \\ f {s)dB^'^'> {s)ds\\ . 2) = WfWh ^'^^ ^11 / 



since we 



have||:7''(/)||(^.) 



I(L2) 

for such /. Since Wiener integrals with respect to standard 



Brownian motion are the elements of the set {J^ /(s) dB{s), f E i^(R)} = {/j^ /(s) dB{s), f e 
it seems natural to give the following definition. 

Definition 5.2. (Wiener integral with respect to mBm) 
For an mBm B^'^\ define the Gaussian space by 



e, :={/ f{s)dBi'^){s), feem 

We call Wiener integral with respect to S'-'*) the elements of Oh- 



Remark 5.26. (i) Obviously, Qh = {ES (/r /(s)il[5e, (s, < ■,ek> 
to definition \4-l\ theorem \5.25\ and the fact that J^{f) = /(s) dB^^\s) on £{ 



f e Sim 
[), we have 



Thanks 



(5.47) 



In other words, the set of Wiener integrals in the sense of definition \4-.l\ and \5.2\ coincide. 

(a) When h is a constant function equal to H we find that Qh ~ Qh = {< -jMnif) > : / <= L'jj{M.)} 

since £{M.) ' = £{M.) ' " = i|^(R). This is exactly what is expected in view of l\4-^ - 

In fact we can be a little more precise in the case of fBm. Let supp{K,) denote the set of measurable 
functions / : M — > R with compact support. 

Proposition 5.27. Let H e (0, 1). Then: 

(i) Let f ■.'R^R he in L]^^{R) n L]j{R). Then S^f{s) dB^"\s) belongs to (L^) if and only if 
J^fis) dB("\s) = j"{f). 
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(ii) i^W n ilf(K) n suppiK.) C {/ : M ^ R : /r/(s) dB("\s) G (L^)}. 

I'm; For fi-almost every f m i,^„^(R) n supp{K) n {/ : M ^ R : /jj/(s) dB^"\s) € (i^)}, / is m 
L2^(R) and verifies /r/(s) dB^"\s) = 

Proof, (i) Let / e L\^^{M)r\B]j{M) and define $/ ;= dB^"\s). By theorem [Q^ J2k=a < 

fjMniek) >< -jCfe > where the equahty holds in (S)* . If we assume that $/ belongs to (L^), then the 
equality is valid in (i^). Besides, since / belongs to L;\,^(R) n L'jj{M.) we have, according to theorem 13. 71 
J" if) .,MH{f) ES< MH(/),efc >i.(R) < .,efe >= < LMniek) >< .e^ > in (L^). 

The converse part is obvious since /^/(.s) dB^"^s) = entails that /,j/(s) dB'-"'>{s) belongs to 

(ii) Since / is in Ll^^{K.) n supp{JC) theorem 15.251 entails that 

j^f{s)dBW^s) = E+^o (/R/(s)M^(efe)(s)ds) < .,e, > = < f^Mnieu) >< ..e^ > in (S)* (one 

only needs to replace il[o,t] by ilsupp(/): where supp{f) denotes the support of / in theorem I5.25p . 
Besides, since / belongs to L|^(R), J^{f) exists and is equal to Yl,k=o < fT^H^ek) >< -^ek > in (L^)- 
{Hi) Fix / in Ll^{M.) n supp{K) n {/ : M ^ R : /r/(s) dB^"\s) £ (i^)} and define h as subset of w 
in n such that (|5.8p is true for all g in {cfc : fc G N}. Clearly i7 belongs to 5. As soon as / is in H, we 
can write /r/(s) dB^"^s) = J2k=o < MH{f),ek >< .,ek > in (L^). This entails that Mnif) belongs 
to L^{R) and then, by bijectivity of Mh, that / belongs to Ljj(R). □ 

Remark 5.28. This proposition shows in particular, that for fi— almost every g in supp{lC): 

5Gi^(R)^ / g{s)dB("\s)£{L^). 

Moreover, in this case, 

J"i9) ^= I 9{s)dBi"\s). 

6 It 6 Formulas 

6.1 Ito Formula for generalized functionals of mBm on an interval [a, 6] with 

< a < 

Let us fix some notations. For a tempered distribution G and a positive integer n, let G^"-* denote 
the n*'^ distributional derivative of G. We also write G" :— G^^\ Hence, by definition, the equality 
< G', If >= - <G,tp' > holds for all tp in =5^(R). For a map t F{t) from [a, h] to J^_p(R) we wih note 

^§^{t) the quantity (F(t))^"\ that is the n*-^ derivative in ^'(R), of the tempered distribution F{t). 
Hence we may consider the map t ^^if) from [a, b] to ,S^'{M). Moreover for any t^ in [a, b], we will 
note ^(to) the quantity lirn -^(*»+'')~-'^(*o) ■^];icn it exists in ^_p(R), for a certain integer p. When it 

exists, ^{to) is a tempered distribution, which is said to be the derivative of the distribution F(t) with 
respect to t at point t = to- In line with section [5.31 we then define, for Iq in [a, b] and a positive integer 
n, the following quantities: 

— {to,B'^{toj) (F(to))<"'(S^to)) and _(t„, i?'^(t„)) (^(to)j (S^to))- 

Theorem 6.1. Leip e N, a andb two reals withQ < a < b, and let F be an element of {[a,b], S'^p(U.)) 
such that both maps ^ and • from [a,b] into J?^-_p(R), are continuous. Then the following equality 
holds in (S) : 



F{b,B^^'\b))- F{a,B^^'\a)) = 
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Remark 6.2. Recall that for all t in [a,b], ■^[Rh{t,t)] = 2 t2/i(t)-i {^^ (^-^ t\nt + h{t)). 

Proof. Wc follow closely [1] p. 97-98 for this proof. First notice that the three integrals on the right 
side of (|6.ip exist since all integrands verify the assumptions of theorem 15.181 According to lemma 12.81 
it is then sufficient to show equality of the S'-transforms of both sides of (j6.ip . It is easy to see that, 
for every rj g the function t i-> 7 (t^'^^^\ . — J* Mh(j^)[ri) [u) du^ is differentiable from (0,6] into 

■5^(K). Using theorem 15.121 and theorem 15. 171 we may write, for t in [0,6]: 



PiF{t,B('^\t))irj) = j^(F{t),^(^t^'^('\.-J^ A4(,)(^) {u)duyj 
'^{t),l[t''^^'\.~ M,(,)(ry) {u)du 



+ 2 t^'^W-i {h'{t) t\nt + h{t)) (pit), |l (^t^Mt)^ _ _ ^^Jj 
- j^MMm (Fit), g (^t^'^W,. - £ M,(,)(,7) (u) du^"^ =: h+h + h- 

Now, h = 'S(^W(5*''Hi)))('7) = s(^^{t,B^^'>{t))^{r]) using theorem [SH (m). Besides, since 7 
fulfills the equality ^ = ^f;^, we get 

I, = t2Mt)-i i int + /j(t)) ^F(t), (^i^Mt)^ _ £ 



Using (m) of theorem 15. 121 

.(t\ ^ I i^^i^) - /"* 

dx 

Finally, we obtain 



dt 



S{F{t,B'^^\tm) = s(^^{t,B^^\t))){^) + s{w^'^\t)o^{t,B^^\t {^) 

[Rhit,t)]s(^^it,B^'^Ht))){rj). 



1 d_ 

2 rfi' 



□ 



In the proof of theorem 16.61 we will need the particular case where the function F(.) is constant, equal 
to a tempered distribution that we denote F. In this case we have the following 

Corollary 6.3. Let < a < b and F be a tempered distribution. Then the following equality holds in 
(ST: 

F(i?W(6))-^^(B('')(a))=^V'(B('')(s))di?(")(,s) + i j' (^^jR,,{s, s)]^ F" (B^'^Hs)) ds. 

Remark 6.4. Of course when the function h is constant on M, we get the R6 formula for fractional 
Brownian motion given in 
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6.2 Ito Formula in (L^) 

In this subsection, wc give two further versions of Ito formula. The first one holds for functions with poly- 
nomial growth but weak differentiability assumptions, whereas the second one deals with C^'^ functions 
with sub-exponential growth. 

6.2.1 Ito Formula for certain generalized functionals of mBm on an interval 

Theorem 16.11 does not extend immediately to the case a = because the generalized functional is not 
defined in this situation, since M/j(il[o.t]) converges to a.s and in when t tends to (see theorem- 

definition 15. 2|) . As in [4] , we now extend the formula to deal with this difficulty. We will need the 
following lemma which is a particular case of lemma 16.81 below. 

Lemma 6.5. Let / : M — > R 6e a continuous function such that there exists a couple (C, A) m M x IR+ 
with \f{y)\ ^ Ce'^y , for all real y. Let g : M. M.*^ be a measurable function such that limg{t) = and 

define Lf on M.*^ xR by Lf{u,x) := J„ f{y)j{u,x — y) dy. Then lim Lf{g{t),x) — f{xQ), for all 

{t,x)^{0+ ,xo) 

real xq . 

Theorem 6.6. Let _F : R — > R 6e continuous at and of polynomial growth. Assume that the first 
distributional derivative of F is of function type {defined at the beginning of section \3.1^ . Then the 
following equality holds in {L"^): 

F(B('')(6))-F(0)=^V'(B('')(s))di?W(s) + i (^^jR„is, s)]^ F'^B'^'^Hs)) ds. (6.2) 
Proof. We follow again closely [3]. 

Step 1: lim F{B^''\t)) = F(0) in (S*). In order to establish this fact, let us use theorem 8.6 of [27]. 

t — ^0-1- 

Since F is of polynomial growth, we may write, thanks to formula (29) of [1], that there exist two reals 
C and M and a positive integer m such that E[FiB'^^\t)f] ^ C^{1 + ^^\tf < M^, for all t in 
[0,b]. Since || : e^ '''^ : ||g — e^'^'lo, Cauchy-Schwarz inequality yields 



|5[i^(i?('')(0)](7?)| = |E[F(i3('^)(t)) : :]| \\F{B^''\tm 



II ■ p<- '^> ■ II 
II ■ -No 



(1 + M) e^^<; for ah t in [0, b] and r] in y{R). (6.3) 

It then just remains to show that lim S[F{B'-'^'> {t))]{ri) = F{0). Thanks to theorem 15.121 and lemma 
16.51 we get 

^hm 5[F(5('^)(t))](,) = ^ln. F(,)-^=i=exp|-^(^* M,,)(.)(u) - .) \ dy 



t-i-0+ 



lim / Fiy) 7 t^'.^*', / M,,^t)iv)iu) du~y) dy 



lim Lf ( / A4(o(^)(^) du ] = F(0). (6.4) 



Step 2: lim^/V'(B('')(t)) dS('*)(t) ^ F' {B^^) [t]) dB'^^\t) in (S)* . 

Define \Hi,H9] [ min h(t), max/i(t)l and let us prove the two following facts 

te[o.b] teloM 

(i) There exists a constant Di which depends only of F such that I |F'(i3('''')(t)||_j^ ^ -Di max{^^,^^} 
for alH in (0,6]. 
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Let us first notice tliat, for all {x, b, t) in M x x (0, b], we have 

exp{-a;V4t2''(*)} g-^'/'' + exp {~x^ /Ab^"^}, 

where e{b) = 1 if 6 ^ 1 and e{b) = if 6 < 1. Note moreover that the function x !—> F'{x) (e^^^/"* + 
e{b) exp {— belongs to L^(M) since F' is of function type and belongs to Since the 

operator A^^ has a norm operator equal to 1/2 (see [57] p. 17) and using the equality |A'/;j(t) (il[o,t])|p = 
t^Kt) ^ we get the following upper bound, valid for all k in N, 



I < F'.i^tMtW > I' = |/m F'{x) 7rV4(fc!)i/2t'=Mt)cxp{-^}efe(x/(y2t''W)) 
7ri/2 sup|e2(^)| (/^ \F'{x)\ exp{-^} dx)' t^^''^^*) fc! 

TT^/^sup I sup|e^(u)| : fc e n} ( / |F'(a;)| (g-^'/^ + cxp {^x"^ / dx) t^fc^t) fcl. 



:-Do 



Using iii) of remark [531 and again the fact that the operator A ^ has a norm operator equal to 1/2 (see 
(2) p. 17 of [57] ) we can write, for all real t in (0, 6], that 



|F'(i?W(i))||'_, 



1 



fe=0 

+ 00 



27rt2'»(t) 



^ (fc!)-^ i-^'^'^W I < > 1^ KA-i)^'^ (( M,(,)(il[o,]) I 



fc=0 



k=0 



= |^-Km,(„(i1,o,,j))| 

\^'^}- (6-5) 



(ii) F'(B('')(t)) dS('')(i) exists in (5)* and is equal to lim / V'(B('') (t)) rfSC') (t) in the sense of 

In order to establish the existence of F' {B^^\t)) dB^^\t) in (5)*, let us use theorem 12.91 From 
theorem-definition l2.21 we know that i^'(i?('')(t))oM^('')(t)) belongs to {S)* for every t in (0, b]. Moreover 
using lemma [2^ we get, for rj e y{W) and t e (0, 6], 



5(F'(S('')(t))oW^W(t))(ry) ^ ||F'(i?W(t))|L, l|W^"'H0ll-2 exp{|r;|^} 

<A"max|^,^| exp{|?7|2}, (6.6) 



= :C(t) 

where wc have defined K := ^ sup ||VF('')(t)|U2- The function t ^ S{F'{B^''^t)) o VF('')(i))(r/) is 

^ te[o,b] 

measurable on [0,b] since S {F' {B'-'^^t)) o W'-''\t)){r]) = S{F' {B'-''\t))){'q)S{W^''^ {t)){r]) using the- 
orems [5321 and [5?TZ1 Moreover, since C belongs to L^{[0,b]), theorem 12.91 applies and shows that 
/q F'(B('')(t)) dB('')(t) is in (S)* . It then just remains to use theorem 8.6 in [57] to show the con- 
vergence, in the sense of (5)*, of /V'(B('')(t)) dB^''\t) to F'(B('') (t)) dB^f^Ht) as a tends to 0+. 
Let {an)ne'M ^ decreasing sequence of real numbers which tends to when n tends to +oo and 
:= F' {B'^^^t)) dB'^^\t) - jl^F'{B'^^\t)) dB'^^\t). For every 77 e J5^(M) and every n e N, 
S{'^n){ii) = /o''il[o,a„](0 SiF^iB^'^^t)) o W^('')(t))(77) dt. Using (HU) and the dominated convergence 
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theorem, it is easy to show that hm S{^n){v) = 0- Hence theorem 8.6 of [37] apphes and shows that 
lim_X!'F'(B('')(t)) = j^F'{B('''\t)) dB^'''\t) in (5)*. 



Step 3: Proof of 

For any real a such that < a < 6, we have, thanks to corohary 

F{B^'\b)) - F{B(^\a)) - j\'{B^'\s)) dB(^\s) ^ \ j' (^^jR„is, s)]J F" (B^'^Hs)) ds. 

Steps 1 and 2 ensure that the left hand side has a limit in (S)* when a tends to 0. Using theorem 
[2Jl it is easy to see that (^[i??,(s, s)]) F"{B''^\s)) ds belongs to (5)*. Hence using the dominated 
convergence theorem and lemma {^IM we deduce that lim (^[i?/j(s, s)]) F" {B^^\s)) ds is equal 

to /q (^[i?/i(.s,.s)]) F"{B^''^s)) ds in (5)*. Since we have proved that, for aU t in [0,6], F{B^''^t)) 
belongs to (i^), the same holds for the right hand side of (|6.2p and then this equality holds also in 
(L2). □ 

Remark 6.7. As in the case of fBm (see J^), the fact that both sides of the equality (jg.^p are in (L^) 
does not imply that every single element of the right hand side is in [L^). This will he true if, for instance, 



F"{B^''^t)) belongs to (i^) and 







f 

'0 


ds 



\F"{B^''\i 



ds < +00. 



6.2.2 Ito Formula in (L^) for C^'^ functions with sub-exponential growth 

Let us begin with the following lemma: 

Lemma 6.8. Let T > and / : [0, T] x M — > K 6e a continuous function such that there exists a couple 
{Ct, At) o/M X M.*^_ such that max \ f{t,y)\ ^ Crc^'^y for all real y. Define a > At, la ■= (0, j^) and 

: M X R_|_ X la ^ M by Jf{x,t,u) := f{t,y)j(u,x — y) dy. Then Jf is well defined and moreover 
lim J(x, t, 7i) = /(O, xq). 

(x,t,u)^{xo,0+ ,0+) 

Proof. This is an immediate consequence of theorems 1 p. 88 and 2 p. 89 of [35] • □ 

Let us now give an Ito formula for functions with subexponential growth. 

Theorem 6.9. Let T > and h : W ^ (0,1) be a function such that h' is bounded on K. Let f be 
a C^''^{[0,T] X R,]R) function. Furthermore, assume that there are constants C ^ and A < 4 „^ax^ t^h{t) 

t£[o,r] 

such that for all (t, x) in [0, T] x R, 



max { \f(t,x)\, 
te[o,T] I '•' ^ ' ^" 



a/ 

dt 



{t,x] 



dl 

dx 



it,x] 



9V 

dx'^ 



{t,x] 



^ Ce 



Xx-' 



(6.7) 



Then, for all t in [0,T], the following equality holds in (L^): 



/(r,B('')(r)) = /(0,0) + 1^ ^{t,B^^\t))dt + ^.{t^B^'^\t))dB^'^\t) 



Proof. Our proof is similar to the one of theorem 5.3 in [5]. Let T > and t G [0,T]. Formula (|6.8 
may be rewritten as 



dl 

dx 



(t,i?W(t)) dB^'^\t) = /(r,i?('')(r)) - /(o,o) - 



dt 
d 
di 



it,B'^''\t)) dt 



[Rh{t,t)]) ^{t,B(>^\t))dt. (6.9) 
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In order to show that all members on the right hand side of (|6.9p are in (L^), let us use theorem l2.10l If 
G belongs to |/, ^Ji^ we may write, thanks to (|6.7p . 



E 



- 4A 1 



dv 



< (1-4A maxi^M*)) 
te[o,T] 



(6.10) 

Since B^'^^t) = S^(t)|ij=;i(t) ^-s, it is easy to see (in view of [S] p. 978) that i?('')(i) is a Gaussian variable 
with mean equal to Mm){v){u) du = g,,(t, /i(t)) and variance equal to t^'*^*) under the probability 
which has been defined in (|2.18p . Hence, for every t in (0,T] and 77 in ^^(M), 



5(G(i,i?W(t))(7?)=EQjG(t,B('')(t))]= / G(t,u + .g,(t,/i(0)) ^{t^^^'\u) du 

Jr 

G U u + y" A'4(t)(?7)(a;) da; j -i= e'^'/^ ^^^^ 



.11) 



Using the theorem of continuity under the integral sign, we see that the functions 1 1— >■ 5[G(i, B'''^\t))]{ri) 
and t ^ S[G{t, B^'^^t)) o W^'''> {t)]{'q) are continuous on [0,T]. Moreover, in view of ^JU^, t ^ 
\\G{t, B'-^^t))\\Q belongs to L^{[0,T]). Furthermore, note that 



mt,t)]) g(t,BW(t)) 



dti^ 2M / i2h{t)--i \h'{t) t\nt + h{t)\ dt<+c 



Thus, theorem 12.101 applies and shows that all members on the right side of (j6.9p are in (L^). 
Let us now show that 1 1-> ^{t,B'^^\t)) oW'^^\t) is (5)*-integrable over [0,r]. 

Reasoning as in the estimate (|6.6[) . we note that there exists an integer q ^ 2 such that f|(i, B'^^\t)) o 
W^''^t) belongs to (5_,) for every t in [0,T]. Moreover, for every rj in S^{M.) and every t in (0,r] we 
have, using lemma [2761 



s(^^it,B('^Ht)) oW^'^Ht)) (ry) 



^ sup llM^^^Hi)!!- 



te[o,T] 



I — 
9x 



(t,i?('')(t))|| exp{h|^}. 



Since t ^ 1 1 (i, B('')(t)) | j^belongs to ^^([O,^]), theorem HH applies on [0,r] to the effect that 

§|(i, B^^\t)) dB''^\t) belongs to {S)* . It then just remains to show the following equality for all t 
in [0, T] and 77 in 



S ^ g(i, B^^\t)) dB^^\tyj (r;) = S (^/(T, i?('')(T)) - /(0, 0) - ^ |[(t, i?(")(0) dt j (ry) 



Using (|6.1ip with G ~ f and applying the theorem of differentiation under the integral sign, we get 



J^[5(/(t,i?W(t))(^)] =4 i [fit,u + g,{t,hit))) jit^'^i^lu)] du 

^{t^'^W,u)^ [fit,u + g^it,h{t)))] du+ [ fit,u + g,it,hit))) ^ [-/{t^'^^'Ku)] du . 



= .Ui{t} 



= :C/2(t) 



Now, 
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Ulit) = ^ 7(t"'(*', ") ^{t, u + h{t))) du + l^ lit"''-'\u) u + g,it, hit))) j^g.it, h{t))du 

= S [^{t, BC^Ht))) (^) + S BC^Ht)))) S {W^'^t)) (^). 

Besides, using the equahty ^ = ■^^^ and an integration by parts, we get 

Hence we obtain, for any e > 0, upon integrating 1 1-^ Ui{t) + U2{t) between e and T, 

Sif{t,B^'^Ht))irj) - 5(/(£,i?(")(£))(,y) (^^(t,B('^Ht))) iv) dt 

+ [ s[^{t,B^'^\t))oW^^Ht)){r,)dt + l£ ^^[t^Ht)^s[^{t,B(^\t))){^)dt. (6.13) 

Let us now show that hrn^ S (/(e, B^'''>{e)) (77) = /(O, 0) ^ S (/(O, B('')(0)) (77). For every e > 0, (UTTTIl 
can be rewritten as 5(/(e, (e))(77) = ^ /(e,y) 7(£2M^), /^(e)) _ y) rfy. 

For a fixed T > 0, let \t,Ct be such that (|6J)) is fulfilled. There exists b > such that e^'"'^"^ 
belongs to la (defined in lemma as soon as < e < 6. Hence we may write, for any e in (0,6), 
S{fie,B^'^^e)){r]) = J/(g„(e, /i(e)), e, ^^''(e)). Since lim e^'^le) and lim gJe,h{e)) are equal to 0, 

lemma 15^ applies with xq = 0, and yields lim 5* (/(e, i?'-'*' (e)) (?/) ~ /(0,0). 

Let us now establish (|6.12p . Thanks to the fact that both the functions t S[G{t, B^'^^ {t))]{r]) and 
t I— >■ S'[G(t, i?'-'*' (t)) o W^'^'^ {t)]{ri) arc continuous on [0,r] and using the dominated convergence we can 
take the limit when e tends to on the right hand side of (|6.13l) and finally get 

S|/(T, B<">[T)) - /(0,0)|(.,) = S (£ ^{t, Bl*>(i)) dt) (,) + S (^ |i {i. S<«(t)) dB>''Ht)] (,) 

□ 

Remark 6.10. We observe that if we take expectations on both sides of Ito's formula ( | 6. 7p , we get 
exactly formula (1) of theorem 2.1 of f23\j . which is a general weak ltd formula for Gaussian processes, 
in the particular case where the Gaussian process is chosen to be an mBm. 

7 Tanaka formula and examples 

In this section we first give a Tanaka formula as a corollary to theorem 16.61 with F : x n- |a; — a|. We 
then consider the case of two particular h functions that give noteworthy results. 

7.1 Tanaka formula 

Theorem 7.1 (Tanaka formula for mBm). Let /i : M -> (0, 1) he of class , a G M and T > 0. The 

following equality holds in {L^): 
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|BW(r)-a| = |a|+ / sign(B^^\t) ~ a) dB^'^\t) + ^[R„it,t)] S{,}{B^'^\t)) dt, (7.1) 



where the function sign is defined on M by sign{x) := i1r;^(x) — i1r_ (x). 
Proof. This is a direct application of theorem 16.61 with F : a: i~> |a: — a|. 



□ 



Remark 7.2. That the previous equality holds true in (L^) does not imply of course that both integrals 
above are in This last result will be established in a forthcoming paper. 

7.2 Ito formula for functions h such that -^[Rh{t,t)] = 

If h verifies ■^[Rh{t,t)] = 0, then the second order term ^^{t^B^'^\t)) disappears in Ito formula. The 
formula is then formally the same one as in ordinary calculus. In this case, (j7.ip reads: 



|B('')(T) - a| = |a| + r signiB^^^t) - a) dB^''\t). 
Jo 



Note that the "local time" part disappears in this equality. 

Solving the differential equation ^[Rh{t,t)] — yields two families of functions, denoted {hi^x)^^^, and 



VA > 0, /ii,A : (-00, -e^) U (e^,+(X)) 

t ^ 



(0,1) 

A 

In 1*1 ' 



and 



VA<0, /i2,A : (-e^0) U (0,e^) ^ (0,1) 

i ^ 4tT- 



In order to obtain an mBm defined on a compact interval, we may choose a compact subset of (— oo, — e^)U 
(e'^, +oo) when A > and a compact subset of (— e"^, 0) U (0, e'*') when A < 0. 

Figures 1 and 2 display examples of mBm with functions hi{t) := defined on [e + 10^'^, 100] and 
h2{t) -.^ 1^ defined on [lO^^, 1/e - lO^^]. 





30 40 SO 60 70 



90 100 



0.05 0.1 



.15 0.2 0.25 0.3 



0.35 0.4 



Figure 1: t ^ B^'^^^t) with hi{t) 
[e + 10^'\1001. 



Figure 2: t B^^^\t) with h2{t) 
[1.10-3,1/6-10-3]. 



Int 



Note moreover that lim hAt) = lim hAt) = for i = 1, 2. 
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7.3 Ito formula for functions h such that = 1 

The situation where = 1 is interesting since then Ito formula is formahy the same as in the 

case of standard Brownian motion. As a consequence, Tanaka formula takes the familiar form : 

|S('*)(T) - a| = |a| + r sign(i3('^)(t) - a) dB^'^\t) + C 6{a}{B^''\t)) dt. 
Jo Jo 

Thus, instead of a "weighted" local time as in (j7.ip . we get here an explicit expression for the local time 
of mBm for a family of h functions that we describe now. 
The solutions of the differential equation are given by 

: (c,+c^)\{-l,0,l} ^ R 

. 1 In(t-c) (7.2) 

^ 2 ln|i| ' 

where c e R. 

Recall that he is required to range in (0,1). Denote, for c € R, /c := { < G (c, +oo)\{ — 1, 0, 1} : < 
hcit) < 1}. For c in (-oo,l/4), let h := ti{c) := l^^/IEM and ts := <2(c) := ^+^^^ . Then 1^ is 
explicitly given as follows: 



Vc e {-oo, -2] 

Vce (-2,-1], 
Vce (-1,0), 
Vce [0,1/4), 
Vc ^ 1/4, 



Ic= (l + C,ii)U(t2,+0o), 
Ic = (tl,l +C) U (t2,+0o), 

/c-(ti,0)U(0,l + c)U(t2,+oo), 
/c = (ii,i2)U(l + c, +oo), 
/c = (1 + c, +oo). 



Figures 3 and 4 display examples of multifractional Brownian motion with regularity functions /i3(t) 
iMi^ defined on [2 + 10-3,5] and hi{t) -.^ 5^^^ defined on ^ i^-s _in-3i 



10' 



-10- 




Figure 3: 1 1^ B'^'^^^t) with hsit) 
[2 + 10-3,5]. 



1 In(t-l) 

2 Int 




Figure 4: i ^ B'^^^\t) with /i4(t) Inrriip on 



ri-V5 



10-3,-10 
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Note that the case c = yields the constant function = 1/2, i.e. standard Brownian motion. Moreover, 
since lim hc{t) = 1/2 for every c, we see that the family of functions behaves, asymptotically, like 



>-f CO 



the constant function equal to 1/2. However this does not mean that there is convergence in law of _B('''=) 
to Brownian motion; in fact one needs to scale _B('*<='. For every t in R*, define the process Xt on R_|_ by 

Xt{u) := . Then {Xt{u)] u e R+} — {B{u); u e R+} where B stiU denotes a Brownian 

motion and C denotes convergence in law. 



37 



8 Conclusion and future work 



In this paper we have used a white noise approach to define a stochastic integral with respect to multifrac- 
tional Brownian motion which generaUzes the one for fBm based on the same approach. This stochastic 
calculus allows to solve some particular stochastic differential equations. We are currently investigating 
several extensions of this work. In order to apply this calculus to financial mathematics or to physics, 
it is necessary to study further the theory of stochastic differential equations driven by a mBm. This is 
the topic of the forthcoming paper [55] . 

The Tanaka formula we have obtained suggests that one can get several integral representations of local 
time with respect to mBm. Finally, since mBm is a Gaussian process, it seems also natural to investigate 
the links between the construction of stochastic integral with respect to mBm that we gave and the one 
provided by Malliavin calculus. An extension in higher dimension is also desirable. 
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Appendix 



A Bochner integral 



All the following notions about the integral in the Bochner sense come from [22] p. 72, 80 and 82 and 
from g3 p.247. 

Definition A.l (Bochner integral |27| p.247). Let I be a subset o/R endowed with the Lebesgue measure. 
One says that <&:/—>■ (S)* is Bochner integrable on I if it satisfies the two following conditions: 

1. $ is weakly measurable on I i.e u H-<: :> is measurable on I for every ip in [S). 

2. There exists p G N such that <I>(u) G {S-p) for almost every u & I and u H- ||$(m)||_p belongs to 
L\I). 

The Bochner-integral of ^ on I is denoted $(s) ds . 
Properties A.l. If ^ : I ^ (S)* is Bochner-integrable on I then 

1. there exists an integer p such that 

' |$(,s)||_pds. (A.l) 

2. $ is also Pettis-integrable on I and both integrals coincide on /. 

Remark A.l. The previous property shows that there is no risk of confusion by using the same notation 
for both the Bochner integral and the Pettis integral. 

Theorem A. 2. Letp G N and {^n)nen ^ sequence of processes from I to (S)* such that $„(w) G (S-p) 
for almost every u ^ I and for every n. Assume moreover that $„ is Bochner-integrable on I for every 
n and that 

hm /||$,„(s)-$„(s)||_^ds = 0. 

(n,Tn)— i.(+oo,+oo) J j 

Then there exists an {S) -process (almost surely {S-p) -valued), denoted $, defined and Bochner-integrable 
on I such that 



J $(s) ds 
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lim / ||$(s)-$„(s)|| ds = 0. (A.2) 

Furthermore, if there exists an (S)* -process, denoted ^, which verifies l\A.2^ . then ^(.s) = $(s) for 
almost every s in I . Finally the following equality holds: 

lim / $„(s)ds= / $(s)ds m{S*). 

B Proof of proposition 14.21 

Let B^'^^ be a normalized mBm on R with covariance function noted R^- It is well known that 
one can define on the linear space span]g{i?/j(t, .) : t G R} an inner product, denoted <,>f; , by 
< Rh{t, .),Rh{s, .) := Rhit,s) (see [35] p.f20 ff.). Define Eh the closure of spanR{i?,,(t, .) : t'e R} 
for the norm || The space {Eh, \ \ Wr^) is called the Cameron-Martin space (or Reproducing Kernel 
Hilbert Space (R.K.H.S.)) associated to the the Gaussian process S^''^ Let £(R) denote the quotient 
space obtained by identifying all functions of f (R) which are equal almost everywhere. On £'(R) x £(R) 
define a bilinear form, noted <, >;i, by < il[o.t], il[o,s] >/j Rh{t,s). Then <,>/j is an inner product 
provided the linear map Kh '■ £(R) — > Eh defined by Kh{'^^[Q,t]) Rh{t,.), t G R, is injective. De- 



fine Ih := vcctM{BW{t) : t e R} the first Wiener chaos of B^^\ It is a well-known property of 



R.K.H.S. that the map Th : (S^, || ||^^ ) ( spa.n^{B'^'^') (t) : i e R} , || 11(^2)), defined for aU real t by 

ThiRhit, ■)) ~ B^'^\t) is an isometry. As a result, Kh is injective if and only if Th o Kh is injective. The 
next proposition states that this is indeed the case for any continuous function h: 

Proposition B.l. Let h be a continuous function defined on R and ranging in (0,1). The family 
(i?'-''-'(t))jg]g, is linearly independent on R, i.e for every positive integer n, (/3i,/32 • • • ,/?«) in R" and 
(ti,t2 - ■ ■ , tn) in (R*)", such that ti ^ tj for i ^ j, the equality 

n 

Y^P,B^^\t,)=0 a.s, (B.l) 

implies /3i = = • • • = = 0. 

The proof of this proposition requires the following lemma, the proof of which is easy and left to the 
reader. 

Lemma B.2. Define, for t e R, the function At : R ^ C by AtiO := «/ C 7^ and At{0) := t. 

Then, for all t, At is C°° on R and verifies, for every n G N, 



a[-\0)='-^. (B.2) 
n + I 

where A[^^ denotes the n*'' derivative of At- 

Proof, of proposition IB.li Let us use a proof by contradiction. By decreasing n if necessary we 
may always assume that {^1,132 ■■■ , /3n) belongs to (R*)". Besides, thanks to and lemma (0> 
equality (|B.1|) also reads < ., J2]=i f^J -^^'i(tj)(il[o.tj]) >= a.s. By taking Fourier transforms, we get 

ft ^^H^M^,]) "=0 Using (EH) this yields 

n 

where we have defined, for j in {1; 2; • • • ; n}, aj := /3j {ch(tj)) ^ ■ By re-arranging if necessary the {ti)i, 

we may assume without loss of generality that h{ti) ^ h{t2) ^ ■ ■ ■ ^ h{tn). Let card{A) denote the 
cardinal of the set A. We distinguish three cases. 
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First case: card({/i(ti); ft-(i2); • • • ;h{tn)}) = 1. 

Since h{ti) = h{t2) = • • • = H, wc get, by multiplying equality (|B.3p by |^|^^^^^ and taking in- 

verse Fourier transform, X)j=i ^3 ^-'^[o.tj] ~ almost everywhere on R. This entails that {ai; a2] ■ ■ ■ ; Ofn} 
and then /32; ■ ■ ■ ; is equal to {0}. 



(B.4) 



Second case: h{ti) > h{t2)- Using that il[o,i](0 ~ ^t{£,)^ (jB.3p reads: 

n 

By lemma IBT2] and taking the limit when ^ tends to in (|B.4[) . we get ai = which constitutes a 
contradiction. 

Third case: h{ti) = h{t2). 

There exists an integer r in {2; 3; • • • ; n — 1}, (fci, ^2, ■ ■ ■ , fc^) in (N*)'^ with 2 ^ ki < k2 < ■ ■ ■ < = n, 
such that 



h{ti) = /i(i2) = ••• = /i(ifej 
h{tk^+i) = h[tk^+2) = • • • = /i(tfc2) 



= : 



(B.5) 



where 1 > Hi > H2 > ■ ■ ■ > Hr > 0. Note that when r = 1 we have card{{h{ti); ■ ■ ■ ;/i(t„)}) = 1 
(treated in the first case) and when r = n we have card({h(ti); • • • ; /i(t„)}) = n and then h{ti) > h{t2) 
(treated in the second case). We hence assume from now that 2 ^ ki ^ 71 — 1 and 2 ^ r ^ n — 1. 
Define the sets /i, /2, • • • , /,■ by h := {1; 2; • • • ; ki}, I2 := {fci + 1; fci + 2; • • • ; ^2}, ■ ■ ■ , Ir ■= {kr-i + 
1; kr-i + 2; • • • ; fc^}. Using lemma [R2l equality ()B.3[) can be rewritten as 



(=1 je/. 



Lemma B.3. For every p in 



E 



0. 



(B.6) 



Let us admit this lemma for the moment. The equalities {Lp ) for p in {1; 2; • • • ; /ci} yield the following 
linear system: 



h t2 



tki-1 tki 
'-fe- 



fci ,fcl ,fcl ,fcl 

''2 ■ ■ ■ ''fci-1 ''fci 





(B.7) 



The determinant of this system is a Vandermonde determinant which is non zero since all the i,; are 
distinct from each other. As the result, the only solution is ai = (12 = • • • Ofci = which constitutes a 
contradiction and proves the proposition. □ 



We now present a sketch of proof of lemma IB. 31 

Proof of lemma lB.31 By multiplying both sides of equality (|B.6|) by and then taking the 

limit when ^ tends to 0, we get, using lemma [B. 2 [ X^jG/i '^J ^3 ~ ^' which is equality (Li). Now, fix p 
in N* . Starting from equality (jB.6p we first 

- multiply both sides of equality (jB.6|) by |^|^^ '^^^ and call (|B.6I his] the resulting equality. For any ^ 
in M*, we then take the p*'' derivative of both sides of equality (|B.6I bis) at point ^. We call {Ei) the 
equality thus obtained. It reads: 



EE ^3m"''"'At,m 



(p) 



0, V^eM*, 



(El) 
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where [5(C)] denotes the p^'^ derivative of the p-times differentiable map ^ i-^ g{C)- 

Now, starting from \Ei\ , we recursively perform the fohowing operations successively for ? = 2, . . . , r: 

- multiply both sides of equality (-Ej-i) by and call bis) the resulting equality. 

- take the p*^ derivative of both sides of equality {Ei^i bis) at every point ^ in R* and call {Ei) the 
resulting equality. 

Equality {Er) then reads: 



(P) 



•lei 



H3-H2+P 



ip) 



(«) 



Lemma IB. 21 yields that lim aj A[^^ (^) = 



-J-O 



p+1 



= 0, v^eM*, (s,) 

We want to let ^ tend to in the previous 



equality. However, for {l,j) in {1; 2; • • • ; r} x //, lim [• • • ]; /(^) = +00. Nevertheless, it is easy to show 

5— >o 



that, for every in {1; 2; • • • ; r} x lim -^^^^ ^^^h^-h, 



exists in C* and is independent of j . Define 



Q := lim 



■j^:;:z7r, l = l,...,r. Denote, for in {1;2; • • • ;r} x Ui.j 



5^0 All' kr--"' 

map on M such that [• • • = q j^l^-— ^' ^^''^(C) (1 + t^ij(C))- Equality (i;^) then reads 



the continuous 



(B. 



Upon multiplying both sides of the previous equality by ^' we get, for ^ in 



Since ffi > iJ; for Hn {2; • • • ; r}, taking the limit when ^ tends to in l|B.9p and using lemma IBT2] vields 
Ci «P {p + 1)~ X^je/i "^J ^j^^ ~ which is nothing but (Lp+i). This ends the proof. □ 

Remark B.4. Another way to establish that Rh{-, ■) defines an inner product on £(R) for H in (0, 1) 
is to use (|3.3p and p.l3p . 
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